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Linear  Algebra  Applied  to  Physics 


Determining  Small  Vibrations  in  Conservative  Elastic  Systems 


Linear  Algebra  Prerequisites:  Being  able  use  eigenvalues  and 
eigenvector  to  diagonalize  a  symmetric  matrix. 

Prerequisite  Knowledge  in  Physics:  None. 

Other  Prerequisite  Knowledge:  A  background  in  solving  basic 
differential  equations  would  be  helpful.  However,  Appendix  A 
contains  this  basic  information. 
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Section  1  Introduction 

In  this  study  we  will  look  at  small  vibrations.  In  particular, 
the  small  vibrations  which  we  will  study  are  in  a  system  with  an 
equilibrium  configuration  which  is  a  position  where  the  system 
remains  at  rest.  An  example  of  a  system  in  its  equilibrium 
configuration  is  the  simple  pendulum  as  seen  in  Figure  1.1.  The 
simple  pendulum  consists  of  a  ball  attached  to  a  taut  wire,  anchored 
above,  which  can  swing  in  the  vertical  plane.  The  weight  of  the 
wire  is  negligible  compared  to  the  weight  of  the  ball. 

'■//////////y 


Figure  1.1 

We  say  that  a  system  has  a  stable  equilibrium  configuration  if 
after  a  small  displacement,  the  system  tends  to  return  to  Its 
equilibrium  configuration.  There  are  different  types  of  equilibrium 
depending  on  the  nature  of  the  system.  We  are  interested  in  the 
type  of  equilibrium  found  in  an  elastic  system.  This  is  a  system 
which  has  the  following  two  characteristics:  1 )  the  system  has  a 
stable  equilibrium  configuration  and  2)  a  small  displacement  from 
equilibrium  creates  forces  which  tend  to  restore  the  system  to  its 


4 


stable  equilibrium  configuration.  A  displacement  from  equilibrium 
is  called  strain  and  the  force  which  restores  the  system  to 
equilibrium  is  called  stress.  Thus,  stress  is  a  function  of  strain. 

The  simple  pendulum  in  Figure  1 . 1  is  also  an  example  of  an  elastic 
system  in  its  stable  equilibrium  configuration. 

The  total  energy  in  an  elastic  system  is  composed  of  two  types 
of  energy,  kinetic  and  potential.  We  will  begin  by  considering  the 
intuitive  definitions  of  these  terms  and  then  discuss  their  formulas. 
Kinetic  energy  is  the  energy  a  body  possesses  because  it  is  in 
motion.  Before  we  can  write  the  formula  for  kinetic  energy,  we 
must  be  able  to  describe  the  system  mathematically.  In  any  system 
there  is  a  minimum  number  of  coordinates  that  are  required  to  fully 
describe  the  configuration  of  the  system.  In  general,  the  number  of 
coordinates  is  equal  to  the  number  of  “particles”  in  the  system  times 
the  dimension  of  the  system.  In  the  case  of  the  simple  pendulum, 
the  ball  is  the  only  particle  in  the  system.  The  dimension  of  the 
system  is  one,  because  the  position  of  the  ball  can  be  described 
using  the  angle  made  by  the  pendulum  compared  to  the  position  of 
the  pendulum  in  its  equilibrium  configuration,  as  seen  in  Figure  1.1. 
Therefore,  the  number  of  coordinates  needed  to  describe  the  simple 
pendulum  is  one.  The  velocity  of  the  system  can  also  be  written  in 
terms  of  the  coordinates  which  describe  the  configuration  of  the 
system.  To  be  able  to  do  this,  we  must  specialize  our  notation.  If  n 
coordinates  (X|,  X2 . are  required  to  describe  the  system,  then 

each  Xj  represents  a  Cartesian  coordinate  of  one  of  the  particles  in 
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the  system.  For  example,  if  we  have  two  particles  moving  in  the  xy- 
plane,  which  has  dimension  two,  we  will  need  four  coordinates  to 
describe  the  system.  The  four  coordinates  (Xj.  X2.  x^,  x^)  represent 

the  Cartesian  coordinates  of  the  particles  in  the  system;  that  is,  Xj 
and  X2  represent  x-  and  y-coordinates  of  the  first  particle,  and  x^ 

and  x^  represent  the  x-  and  y-coordinates  of  the  second  particle. 

From  this  we  see  that  the  velocity  of  the  system  can  be  expressed  in 
terms  of  the  velocity  of  each  coordinate.  The  velocity  vector  for  a 
system  with  n  coordinates  can  be  written  in  terms  of  its  velocity 
components 


<  dt  ’  dt  .  dt  J 


The  kinetic  energy  of  the  system  is  equal  to  the  sum  of  one  half  the 
square  of  each  velocity  component  times  the  mass  of  the  particle 
which  the  coordinate  describes.  If  we  let  T  represent  kinetic 
energy  and  mj  the  mass  of  the  particle  which  is  described  using  the 

Cartesian  coordinate  Xj,  then  our  formula  becomes 
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To  have  a  conservative  system,  there  must  exist  a  function 
whose  partial  derivative  with  respect  to  any  coordinate,  say  Xj.  is 

equal  to  the  negative  value  of  the  force  in  the  direction  represented 
by  that  coordinate.  This  function  is  called  the  potential  energy 
function.  We  can  describe  the  relationship  between  this  function 
and  the  forces  in  the  system  by  the  equation 

■5^(potential  energy  function)  =  - (force  in  the  x. direction) 
dX .  I 

1 

From  now  on.  we  will  assume  that  we  are  always  in  a  conservative 
system.  In  addition,  if  the  potential  energy  function  is  not  time 

dependent  ^  ^(potential  energy  function)  =  0  ).  then  in  our 

conservative  elastic  system  the  total  energy  of  the  system  is 
constant  and  is  the  sum  of  the  kinetic  and  potential  energies.  Also, 
when  the  strain  of  the  system  is  zero  (the  system  is  in  its 
equilibrium  configuration,  so  Xj=0  for  all  i).  then  the  partial 

derivative  of  the  potential  energy  function  with  respect  to  any 
variable  must  equal  zero.  This  statement  can  be  interpreted  in  the 
following  two  ways:  I )  in  the  equilibrium  configuration  the  potential 
energy  function  is  at  a  minimum  and  2)  the  restoring  forces  are 
equal  in  magnitude  and  of  opposite  sign  to  the  forces  that  created 
the  displacement.  This  statement  also  tells  us  that  the  potential 
energy  function  can  not  contain  linear  terms  which  have  nonzero 
constant  coefficients  in  any  of  the  X|.  To  see  why  this  is  true,  let  us 
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assume  the  potential  energy  function  contains  a  nonzero  linear  term 
cXj  (c  is  a  nonzero  constant).  Then  take  the. partial  derivative  of  it 

with  respect  to  Xj.  Setting  Xj  equal  to  zero,  we  find  the  nonzero 

constant  c  is  equal  to  zero,  which  is  a  contradiction.  Therefore,  we 
conclude  that  c  must  be  zero  and  the  potential  energy  function  does 
not  contain  a  nonzero  linear  term  cXj.  Also,  it  does  not  matter  if  the 

potential  energy  function  has  a  constant  term  or  not,  because  when 
we  differentiate  the  function  with  respect  to  any  X|  (i=  1 , 2 . n),  the 

constant  becomes  zero.  Thus,  if  we  write  the  potential  energy 
function  in  its  Taylor  series  expansion,  the  non-constant  part  starts 
with  quadratic  and  terms  of  higher  powers  (which  may  also  contain 
a  constant  term).  When  we  differentiate  the  potential  energy 
function  with  respect  to  Xj  (for  i=l,  2,  ...  .  n),  we  obtain  a  linear 

combination  of  the  variables  x  j ,  X2 . x^^  plus  higher  order  or 

mixed  terms  (for  example  XjX^  or  XjX2X^).  If  we  ignore  the  higher 

order  terms,  then  the  linear  part  which  remains  gives  us  the  specific 
relationship  of  stress  to  strain,  which  is  known  as  Hooke's  Law.  In 
general,  Hooke's  Law  states  that  “stress  is  a  linear  transformation 
operating  on  strain.”  Intuitively,  we  would  say,  restoring  forces  are 
linearly  proportional  to  the  displacement  of  the  mass  from 
equilibrium.  If  the  non-constant  potential  energy  function  starts 
with  a  power  greater  than  two,  it  is  possible  to  use  an  approximation 
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to  find  the  relationship  between  stress  and  strain  in  which  the  still 
higher  power  terms  in  the  partial  derivative  of  the  potential  energy 
function  have  been  ignored.  However,  this  is  no  longer  a  linear 
function. 

We  will  consider  two  approaches  to  the  formulation  of  a 
differential  equation  whizh  models  a  system.  The  first  approach  is 
developed  using  Newton  s  second  and  third  laws  of  motion,  which 
are  stated  below  for  convenience. 

Second  Law  The  mass  of  the  body  times  the  acceleration  of 
the  body  is  equal  to  the  force  acting  on  the 
body. 

Third  Law  For  every  action  there  is  an  opposite  and  equal 
reaction. 

From  these  laws  we  derive  the  differential  equation  which  models  a 
conservative  elastic  system 

2 

d  (displacment) 

mass  X - - =  restoring  force 

dt"^ 

or 


9 


In  the  second  approach,  instead  of  using  the  direct  application 
of  Newton's  laws,  we  will  consider  a  method  developed  by  Joseph 
Lagrange,  a  French  mathematician.  This  very  elegant  and 
sophisticated  method  can  be  applied  to  systems  which  are  more 
general  than  the  ones  we  are  considering  here.  Since  our  system  is 
conservative  and  elastic,  the  energy  is  constant  and  equal  to  the  sum 
of  the  kinetic  and  potential  energies.  If  we  let  V  represent  potential 
energy  and  E  represent  the  total  energy  in  an  n  coordinate  system, 
then  we  have 


E  =  T^V  =  iS  m,[^J 

I  =  I 


+  V 


To  make  it  easier  to  express  this  differential  equation,  we  will 
introduce  a  type  of  notation  you  may  not  have  used  before.  The 
dx 

derivative  —  will  be  written  as  x.  where  the  single  dot  above  the 

variable  x  indicates  that  one  derivative  of  x,  with  respect  to  time, 
has  been  taken.  This  idea  can  be  extended  so  that  x  indicates  that 
two  derivatives  of  x  with  respect  to  time  have  been  taken.  This 
notation  is  used  to  rewrite  the  equation  for  total  energy. 

E=T  + V=y  +V 


We  wish  to  derive  the  equations  of  motion  which  can  be  used  to 
model  conservative  elastic  systems.  Our  first  step  is  to  find  the 
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partial  derivative  of  Equation  (1.2)  with  respect  to  each  of  the 
coordinates.  Since  the  procedure  is  the  same  when  taking  the 
partial  derivative  with  respect  to  each  coordinate,  we  will  only  find 
the  partial  derivative  of  the  function  for  total  energy  (a  constant) 
with  respect  to  the  Xj  th  coordinate.  We  obtain 


0  = 


ar  .  dv 
ax.  ax. 

I  I 


or 

(1.3) 


av  _  dT 
ax.  ~  ax. 

I  > 


Since  the  mass  of  each  particle  is  known,  the  partial  derivative  of 
the  kinetic  energy  with  respect  to  the  coordinate  Xj  is 


ax 

ax. 


-m. 


3x. 
2  X  — - 


=  m 


ax. 

_ I 

at 


ax. 

_ [ 

ax. 


ax. 

=  ni.  =  m  .x. 

I  dt  II 


Substituting  this  into  Equation  (1.3),  we  obtain  the  restoring  force  of 
the  X|  coordinate. 


av 

ax. 


=  m 


I 


(1.4) 


1 


Since  the  kinetic  energy  is  expressed  in  terms  of  x. ,  we  can  find 

3T 
dx.  . 


I 


We  now  differentiate  this  equation  with  respect  to  time  to  obtain 


We  see  that  the  right  side  of  this  equation  is  the  same  as  the  right 
side  of  Equation  ( 1 .4).  Equating  the  two.  we  obtain  the  equation  of 
motion  for  the  Xj  th  coordinate. 

dt  3x.  9x. 

L  I  J  I 

Therefore,  the  equations  of  motion  which  model  our  conservative 
elastic  system  with  n  coordinates  are 


d 

■  ax 

1_  av 

d 

■  ax  ■ 

_  av 

d 

■  ax  ■ 

av 

dt 

.3^1. 

ax,' 

dt 

.^^2. 

3x2' 

■ '  dt 

ax„ 

We  will  model  the  simple  pendulum  of  Figure  1.1  using  both 
the  method  which  applies  Newton's  laws  directly  and  the  equations 
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of  motion  formulated  by  Lagrange.  Since  we  are  only  interested  in 
small  vibrations  of  a  system,  let  us  discuss  the  conditions  under 
which  the  vibrations  in  the  pendulum  system  remain  small.  In  a 
conservative  elastic  system  the  total  energy  of  the  system  is 
constant  and  is  the  sum  of  the  kinetic  and  potential  energies.  The 
potential  energy  of  the  pendulum  system  is  determined  by  the 
displacement  of  the  ball  from  its  equilibrium  position.  Imagine  the 
pendulum  in  Figure  1.1  being  placed  very  close  to  its  equilibrium 
position  and  released  as  in  Figure  1.2. 


y/////////y 


Since  the  potential  energy  is  small  to  start  with  (the  displacement 
from  equilibrium  is  small)  and  we  are  in  a  conservative  system,  we 
know  that  it  will  remain  small.  Because  the  displacement  stays 
small,  the  angle  0  will  always  be  small.  Thus,  the  vibrations  of  this 
system  can  only  be  small  vibrations. 

Let  us  model  the  simple  pendulum  system  using  the  method 
which  applies  Newton’s  laws  directly.  To  keep  this  example  simple 
we  will  only  consider  the  positive  region  which  is  to  the  right  of  the 
equilibrium  configuration  in  Figure  1.3(a).  Let  L  represent  the 
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length  of  the  pendulum,  m  be  the  mass  of  the  ball  at  the  end  of  the 
pendulum,  0  the  angle  the  pendulum  makes  with  respect  to  the 
equilibrium  configuration,  and  s  represent  the  length  of  the  arc  the 
ball  travels.  The  force  pulling  the  ball  down  is  the  mass  of  the  ball 
times  the  acceleration  due  to  gravity  g. 


Figure  1 .3 


We  will  apply  Newton's  laws  of  motion  to  model  the  pendulum 
system  by  using  Equation  (1.1 ).  Thus,  we  need  to  determine  the 
stress  and  the  strain  of  the  system.  Since  stress  or  restoring  force  is 
the  force  trying  to  return  the  ball  to  its  equilibrium  configuration, 
we  must  resolve  the  force  on  the  ball  (mg)  into  its  component  forces. 
Figure  1.3(b)  shows  the  restoring  force  is  the  component  of  force  on 

the  ball  along  the  arc  length.  Since  sin9  = 

**  ®  hypotenuse 

magnitude  of  the  restoring  force  is  mgsin0  as  seen  in  Figure  1.3(c). 

We  will  need  a  minus  sign  to  indicate  that  the  restoring  force  is 
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opposite  in  direction  to  the  force  which  originally  moved  the  ball 
from  its  equilibrium  configuration.  Thus,  the  restoring  force  or 
stress  is  equal  to  -  mgsin0.  The  strain  is  the  displacement  of  the 
ball  from  the  equilibrium  position.  This  distance  is  the  arc  length  s, 
which  can  also  be  described  using  the  equation  s-L0.  Substituting 
these  values  for  stress  and  strain  into  equation  (1.1)  gives 


m 


=  -  mg  sin0 


Taking  the  second  derivative  of  L0  with  respect  to  time,  this 
equation  becomes 


m  L0  =  -  mg  sin0 

Simplifying  and  moving  all  terms  to  the  left  side  of  the  equation,  we 
get  the  second  order  differential  equation  that  models  our 
conservative  elastic  system. 

0  +  -^sin0  =  O  . 

We  now  model  the  simple  pendulum  system  using  the 
equations  of  motion  formulated  by  Lagrange.  Since  0  is  the  only 
coordinate  needed  to  describe  the  system,  we  will  only  need  to  use 
one  of  the  equations  of  motion  found  in  Equation  (1.5). 


I 
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AriIl=_iV 
dltaej  30 

Thus,  we  need  to  find  both  the  kinetic  energy  and  potential  energy 
of  the  system.  The  kinetic  energy  is  one  half  the  mass  of  the  ball 
times  the  square  of  the  velocity.  The  velocity  is  the  first  derivative 
of  the  distance  with  respect  to  time. 

velocity  =-^[ distance ]  =  -^  L0  =L^-L0 

Therefore,  the  equation  describing  the  kinetic  energy  becomes 

T  =  ^m[L0]^=YmL^O^ 


Since  potential  energy  is  the  energy  needed  to  restore  the  system  to 
equilibrium,  it  is  equal  to  weight  of  the  ball  (mass  of  the  ball  times 
the  acceleration  due  to  gravity  g)  times  the  height  of  the  ball  above 
the  reference  point.  Since  the  ball  is  below  the  reference  point,  V  is 


negative.  Using  the  fact  that  cos0  = 


adjacent 
hypotenuse  ’ 


we  determine  the 


distance  of  the  ball  below  the  reference  point  to  be  Lcos0  as  seen  in 
Figure  1.4. 


16 


Figure  1.4 


Thus,  the  potential  energy  is 

V  =  -  mgL  COS0 . 


Now,  we  substitute  the  appropriate  partial  derivatives  of  the  kinetic 
and  potential  energy  equations  into  the  equation  of  motion.  First, 
we  find  the  left  side  of  the  equation  of  motion  by  differentiating  the 
kinetic  energy  equation  with  respect  to  0. 

^  =  4-mL^  (20)  =  mL^0 

O0  ^ 


Then  differentiating  with  respect  to  time,  we  obtain 


_d 

dt 


2  - 
m  L  0 


Second,  the  right  side  of  the  equation  of  motion  is 

dV  r  .  3(cos0)  1  r  .  /  .  ^ 

-*^=-|-mgL  — —  J  =  -  [-  m  g  L  (-sin0)  ]  =  -  m  g  Lsin0 
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Equating  the  two  sides,  the  equation  of  motion  becomes 

2  “ 

m  L  0  =-  m  g  Lsin0 

Simplifying  and  moving  all  terms  to  the  left  side  of  the  equation,  we 
obtain  the  second  order  differential  equation  that  models  our 
conservative  elastic  system.  As  expected,  this  is  the  same  equation 
which  we  found  by  applying  Newton  s  laws. 

(1.6)  0  +  |.sin0  =  O 

Let  us  pause  for  a  moment  and  discuss  the  relationship 
between  the  potential  energy  function  and  the  component  of  force 
tangent  to  the  path  the  ball  travels  (that  is,  in  the  direction  of  arc 
length).  Recall  that  in  a  conservative  system,  the  partial  derivative 
of  the  potential  energy  function  with  respect  to  any  direction,  gives 

the  negative  of  the  force  in  that  direction.  That  is,  •^  =  -  where 

ds  s 

Fg  is  the  force  in  the  direction  of  the  arc  length  s.  First,  we  need  to 

write  the  potential  energy  function  in  terms  of  arc  length  s.  We  will 
use  the  fact  that  s=L0. 

V  =  - mgL  cos^-^^ 

We  continue  by  differentiating  this  with  respect  to  s  to  obtain 
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m  g  sin9 


Thus,  F^  =  -  -  mg  sinS  is  the  restoring  force,  since  the  simple 

pendulum  is  described  using  only  one  coordinate.  (Recall  the 
discussion  following  Figure  1.3.) 

So  far  we  have  found  the  second  order  differential  equation 
which  models  the  simple  pendulum  system  using  two  different 
methods.  Now,  we  are  ready  to  consider  how  Equation  ( 1 .6)  can  be 
solved.  Since  this  equation  involves  sin0,  we  know  it  is  a  nonlinear 
differential  equation  (See  Appendix  A  for  definition.)  One 
technique  used  to  find  the  exact  solution  (if  that  is  possible)  of  a 
second  order  nonlinear  differential  equation  is  to  first  reduce  it  to  a 
first  order  differential  equation.  Recall  that,  the  equation  of  motion 
was  derived  from  Equation  ( 1 .7).  Since  Equation  ( 1 .6)  was  found  by 
using  the  equation  of  motion,  we  can  use  Equation  ( 1.7)  as  our  first 
order  differential  equation. 


(1.7)  E=  T  + V  =  -mgLcosO 

2 

We  begin  by.solving  for  the  squared  derivative  in  Equation 


(1.7). 
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(1.8) 


2(E  +  mgL  COS0) 


mL 


Let  us  pause  for  a  moment  to  assure  ourselves  that  we  could 
legitimately  use  Equation  ( 1 .7)  as  our  first  order  differential 
equation.  We  will  differentiate  Equation  ( 1 .8)  with  respect  to  time 
and  see  that  we  get  Equation  (1.6). 


2  00  = 


m  L 


0  +  mgL(-0sin0) 


When  simplified,  we  see  this  is  Equation  ( 1 .6). 

0  +  -^sin0  =  O 

Taking  the  square  root  of  both  sides  of  Equation  (1.8)  and  recalling 
that  we  are  only  considering  the  positive  region  to  the  right  of  the 
equilibrium  configuration,  we  obtain 

d0  _  /2(E-t-  m  g  Lcos0) 

dt  “  /  2 

-v  mL 


One  way  to  solve  the  differential  equation  above  for  t.  is  to  isolate  dt 
•  • 

on  one  side  of  the  equation  and  d0  on  the  other.  This  technique  is 
sometimes  called  separation  of  variables. 
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dt  = 


mL 


2(E  +  m  gLcosB) 


d6 


Integrating  both  sides  of  the  above  equation,  we  obtain 


t 


m  L 


2(E  +  mgLcose) 


dO 


This  is  an  elliptic  integral  which  can  not  be  expressed  in  terms  of 
elementary  functions.  Thus,  to  get  any  information  about  the 
solution  to  Equation  ( 1 .6),  we  must  resort  to  numerical 
approximations  or  use  the  fact  that  we  are  dealing  with  a  system 
involving  small  vibrations.  We  also  note  that  the  question  of  finding 
the  inverse  function  0  =  0(t)  of  the  function  above  is.  at  best,  a 
numerical  approximation  problem  and  is  not  even  useful  in 
predicting  values  of  0  at  a  given  time  t.  since  we  are  dealing  with 
small  vibrations.  As  a  point  of  interest,  if  we  were  not  considering 
small  vibrations,  then  the  function  t=  t(0)  and  its  inverse  function 
0  =  0(t)  would  be  the  only  tools  with  which  we  could  obtain 
information  about  the  system. 

Using  the  fact  that  we  are  dealing  only  with  small  vibrations, 
we  consider  the  factor  sin0,  which  makes  Equation  (1.6)  a  nonlinear 
differential  equation.  We  can  write  sin0  as  a  Taylor  series 
expanded  about  zero. 


OO 


(-  1) 


sin0=  y  ... 

„To  (2n  +  D! 


-  2n  +  I  „ 

0  =0- 


3! 


5! 
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Since  we  are  considering  only  small  values  of  0,  the  terms  in  the 
expansion  above  which  contain  powers  of  0  are,  in  practice,  ignored 
(a  very  small  number  raised  to  a  power  greater  than  one  becomes 
even  smaller).  Any  time  terms  are  ignored  we  expect  a  certain 
amount  of  error.  To  determine  the  exact  amount  of  error  would 
require  the  same  type  of  calculation  that  it  would  take  to  solve  the 
original  equation.  However,  since  we  are  considering  only  small 
vibrations,  we  are  assured  the  amount  of  error  will  not  affect  the 
resulting  solution.  Thus,  using  the  Taylor  series  expansion  for  sin  0 
we  see  that  sin0  can  be  replaced  by  0.  for  small  values  of  0.  This 
substitution  is  only  valid  when  we  are  dealing  with  small  vibrations. 
Using  this  substitution,  the  second  order  differential  equation 
(Equation  (1.6)  )  becomes 


0  +  |-  0  =  0 


This  is  a  second  order  linear  differential  equation  whose  solution  is 
found  using  basic  techniques  from  differential  equations.  (Basic 
solution  techniques  are  found  in  Appendix  A.)  VVe  obtain 

e  =  C|Cos(yf  t)*C2Sin(yf  t) 

In  summary,  we  have  investigated  two  different  ways  to 
model  a  conservative  elastic  system.  One  method  applies  Newton  s 
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second  and  third  laws  directly  to  the  system  to  create  the 
differential  equation.  The  other  method  uses  a  technique  developed 
by  Lagrange,  which  was  much  easier  to  generalize  and  could  be 
applied  to  many  different  types  of  systems.  The  derived  equations 
of  motion,  greatly  simplify  the  amount  of  work  necessary  to  model  a 
conservative  elastic  system.  From  these  techniques,  we  found  the 
second  order  nonlinear  differential  equation  that  models  the  simple 
pendulum.  Since  we  considered  only  small  vibrations,  we  found 
that  the  equation  could  be  represented  by  a  second  order  linear 
differential  equation  which  has  an  elementary  solution. 


Section  2  Linear  Spring-Mass  Systems 

In  the  last  section  we  considered  a  system  which  only  needed 
one  coordinate  to  completely  describe  the  system.  We  now  look  at 
higher  dimensional  systems,  such  as  spring-mass  systems  in  which 
more  than  one  coordinate  is  required  to  specify  the  state  of  the 
system.  A  spring-mass  system  is  a  conservative  elastic  system  with 
a  stable  equilibrium  position  occurring  when  all  of  the  coordinates 
are  set  equal  to  zero.  To  become  familiar  with  spring-mass  systems, 
we  will  first  consider  the  one  dimensional  case.  Figure  2.1  shows 
the  system  in  its  equilibrium  configuration  (the  spring  is  not  being 
stretched  or  compressed)  where  m  is  the  mass  of  the  block  and  L  is 
the  natural  length  of  the  spring.  We  are  considering  the  spring- 
mass  system  moving  along  a  horizontal  track  rather  than  hanging 


vertically  so  that  we  do  not  have  the  added  complication  of 
describing  how  gravity  affects  the  system. 
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I - L - 1 

Figure  2.1 

To  determine  the  number  of  coordinates  we  need  to  describe  this 
system,  recall  the  formula  given  in  Section  I.  (The  number  of 
coordinates  =  (number  of  particles)  times  (dimension  of  the 
system).  )  The  only  particle  in  the  system  is  the  block  and  since  the 
block  is  moving  along  a  horizontal  track,  the  dimension  of  the 
system  is  one.  Thus,  we  need  only  one  coordinate  x,  to  describe  the 
system.  Imagine  the  block  being  moved  to  the  right  causing  the 
spring  to  be  stretched  x  units.  This  is  shown  in 
Figure  2.2. 


L - 1 - x-H 

Figure  2.2 


To  describe  the  energy  of  this  system  we  again  need  to  find  both 
the  kinetic  energy  and  potential  energy.  The  kinetic  energy  is  one 
half  the  mass  of  the  block  times  the  velocity  squared.  The  equation 
describing  the  kinetic  energy  is 
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T  = 


The  energy  stored  in  the  spring  or  the  potential  energy  of  the  spring 
is  one  half  the  spring  constant  times  the  square  of  the  distance  that 
the  spring  is  stretched.  From  the  laws  of  physics  we  know  that  the 
external  force  acting  on  the  spring  is  proportional  to  the  increase  in 
length  of  the  spring.  We  call  the  constant  of  proportionality  that 
allows  us  to  write  this  relationship  as  an  equation,  the  spring 
constant  or  the  stiffness  of  the  spring  and  each  spring  has  its  own 
specific  spring  constant.  If  we  let  k  represent  the  spring  constant 
and  X  the  displacement  of  the  spring  from  equilibrium,  then  the 
equation  for  potential  energy  is 

V  =  ikx^ 

Since  x  is  the  only  coordinate  needed  to  describe  the  system,  we  will 
only  need  to  use  one  of  the  equations  of  motion  found  in  Equation 
(1.5). 


dt  LdxJ 

To  find  the  left  side,  we  first  differentiate  the  kinetic  energy 
function  with  respect  to  x. 


X 
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Now  differentiate  this  equation  with  respect  to  time. 


To  find  the  right  side  of  the  equation  of  motion,  we  differentiate  the 
potential  energy  function  with  respect  to  x. 


Equating  these  two,  the  equation  of  motion  becomes 

mxs- k  X. 


Simplifying  and  rearranging  terms,  the  differential  equation  which 
models  this  system  is 


m  *  =  “, 

This  is  a  second  order  linear  differential  equation  whose  solution  is 
found  using  basic  techniques  from  differential  equations.  (See 
Appendix  A.)  Note,  the  similarity  between  this  differential  equation 
and  the  one  that  models  the  simple  pendulum. 


X'  =1  “S  {yfk  ‘  )+  (7^  ‘  ) 
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To  be  able  to  model  higher  dimensional  spring-mass  systems, 
we  need  to  study  the  theory  which  describes  the  energy  of  the 
system  in  general  terms.  In  the  one  dimensional  spring-mass 
system  there  was  only  one  coordinate  which  we  labeled  as  x  and  it 
was  expressed  in  terms  of  time.  The  kinetic  energy  of  the  system 
was  described  using  the  first  derivative  of  this  coordinate  with 
respect  to  time,  while  the  potential  energy  was  expressed  in  terms 
of  the  coordinate.  If  we  are  working  with  a  higher  dimensional 
system  which  has  n  coordinates,  say  Xj.  X2 . x^^,  then  the  kinetic 

energy  will  be  described  using  the  first  derivative  with  respect  to 
time  of  each  of  the  coordinates  and  the  potential  energy  will  be 
expressed  in  terms  of  these  n  coordinates.  In  general  we  have 

j  n  rdx."|^  j  n  2 

Kinetic  Energy  ^“7  2)  l“dt^J  “7  S 

i  =  I  i  =  1 

s 

Potential  Energy  ^  =  X  ^  ■  where  V.  is  the  potential 

i  =  I  '  ‘ 

energy  of  each  spring  and 

s  is  the  number  of  springs. 

Since  stable  equilibrium  occurs  when  xi=X2=  ...  =Xn=0,  we  may 
assume  that  the  energy  of  the  system  is  at  a  minimum  in  stable 
equilibrium.  This  means  the  derivative  with  respect  to  any  variable 
must  be  zero  when  that  variable  equals  zero.  Thus,  if  we  have  a 
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function  which  we  wish  to  expand  using  its  Taylor's  series 
expansion,  as  we  did  with  sin0  in  Section  1.  the  expansion  can  not 
have  a  nonzero  linear  term.  For  if  it  did  and  we  took  the  derivative 
of  it,  we  would  end  up  with  a  nonzero  constant.  Subsequently, 
when  all  variables  are  set  equal  to  zero,  the  constant  would  remain, 
indicating  that  we  do  not  have  stable  equilibrium,  a  contradiction. 
Therefore,  the  Taylor  series  expansion  for  the  potential  energy  does 
not  have  linear  terms.  However,  this  expansion  may  have  constant 
terms. 

Let  us  return  to  the  spring-mass  system.  Figure  2.3  shows  a 
system  in  equilibrium  with  two  blocks  having  the  same  mass  m  and 
three  springs  possessing  the  same  length  and  spring  constant.  To 
determine  the  number  of  coordinates  needed  to  describe  this 
system,  we  need  to  recall  the  formula  given  in  Section  I.  (The 
number  of  coordinates  =  (number  of  particles)  times  (dimension  of 
system).  )  The  two  particles  in  the  system  are  the  two  blocks  and 
since  both  blocks  are  moving  along  a  horizontal  track,  the  dimension 
of  the  system  is  one.  Thus,  we  will  need  two  coordinates,  Xj  and  X2, 

to  describe  the  system. 


■■i 

1  T 

1 

T 

1 

T 

Bj 

'  ^1 

1 

^2 

1 

^3 

m 

Figure  2.3 
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Imagine  the  two  masses  are  moved  to  the  right  causing  the  first  two 
springs  to  stretch  by  different  amounts  and  causing  the  third  spring 
to  be  compressed.  This  is  depicted  in  Figure  2.4. 


I — Xj — I 


Figure  2.4 


Now  we  determine  the  second  order  differential  equation  that 
models  this  system.  Thus,  we  need  to  find  the  kinetic  energy  and 
the  potential  energy  of  the  system.  The  equation  below  describes 
the  kinetic  energy  of  the  system  shown  in  Figure  2.4. 


The  potential  energy  of  the  system  is  the  sum  of  the  potential 
energies  of  each  spring.  Spring  1  is  stretched  from  its  equilibrium 
position  by  the  amount  Xj,  so  the  potential  energy  for  spring  1  is 


V|  =s-^k  X  j .  Spring  2  is  stretched  from  its  equilibrium  position  by 
the  amount  X2-X|,  so  that  is  the  potential  energy 


for  spring  2.  Spring  3  is  compressed  from  its  equilibrium  position 
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by  the  amount  X2.  Thus,  the  potential  energy  for  spring  3  is 
^kx  2-  Therefore,  the  potential  energy  of  the  system  is 


v=  i  V 

i  =  1 


.  =  ik[x^^.(x2-x,)  +x'J  =  ik[2x^- 


2X|X2+2X2_ 


Since  this  system  is  described  using  two  variables,  Xj  and  X2.  our 
two  equations  of  motion  are 


First,  we  determine  the  left  side  of  each  equation  of  motion  by 
differentiating  the  kinetic  energy  with  respect  to  Xj  and 

-|I-=im[2x,  +  0]=mx, 

|i  =  lm[0  +  2x2]  =  mx2 

Then,  differentiating  each  of  these  equations  with  respect  to  time, 
we  have 
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dt 


dt 


3T 


ax 


u 


aT 

ax. 


=  •^["’*2]=  "’*2 


To  determine  the  right  side  of  each  equation  of  motion,  we 
differentiate  the  potential  energy  with  respect  to  Xj  and  X2. 


-|^  =  4k[0-2x,+  4xJ=k[x,-2xJ 
Substituting  this  information  into  the  equations  of  motion,  we  obtain 


dt 


ar 

ax 


av 

ax. 


which  becomes 


mX|  =  k[-2X|  +  X2]  or  = -j^ [- 2x ^  +  x J 


and 


dt 


ar  1  av 


ar  1_. 

ax^  ax. 


which  becomes 


X2  =  -^[x,-2x2] 


mx2=k[X|-2x  J 
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We  now  have  a  system  of  second  order  differential  equations  which 
can  be  written  as  the  following  matrix  equation,  where  A  is  a 
symmetric  matrix. 


Since  A  is  a  symmetric  matrix,  all  of  its  eigenvalues  are  real  and  A  is 
diagonalizable.  We  begin  the  determination  of  the  eigenvalues  of 


the  matrix  A  by 

det(A->.I)  =  det  ~  j  =  (- 2 1  =X^+4X  +  3  =  (X+ 3)(X+  l). 


If  we  set  det(A-XI)  equal  to  zero  and  Solve  for  X,  we  find  the 
eigenvalues  are  X=-3  and  X=-l.  Since  A  is  diagonalizable  there  exists 
an  invertible  (orthogonal)  matrix  P  such  that  P  *AP  =  D.  The  matrix 

D  is  the  diagonal  matrix  whose  entries  along  the  main  diagonal 
consist  of  the  eigenvalues  of  A  and  the  columns  of  P  are 
corresponding  eigenvectors  associated  with  these  eigenvalues.  To 
find  P  we  need  to  find  an  eigenvector  associated  with  X=-3  and  one 
associated  with  X=- 1 .  For  X=-3  we  have  the  following  matrix 
equation. 
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!f  we  let  X  j  -  1 ,  then  X2  -  - 1 .  and  it  follows  that  an  eigenvector 

associated  with  the  eigenvalue  X=-3  is  ^  _  J  j.  To  find  an 

eigenvector  associated  with  X=-I,  we  use  the  following  matrix 
equation. 


f  -2+  1 

»  1 

X 

j-r  n 

1  1 

-2+  1  j 

^  2  J 

'1  1  -iJ 

2  J 

U  J 

so  x^ 


If  we  let  x  j  =  1 ,  then  X2  =  1 .  Thus,  an  eigenvector  associated  with 


the  eigenvalue  A,=- 1 


■■{;) 


Therefore,  these  two  eigenvectors  are 


the  columns  of  the  invertible  matrix  P 


■(.; :) 


Recall  that  our  goal  is  to  solve  the  second  order  differential 

equation  X  =  A  X  .  If  we  multiply  both  sides  of  this  equation  by 
P  ^  and  use  the  identity  P  P  ^  =  ^2'  2x2  identity  matrix,  we 

obtain 


(2.1)  P 


‘  x  =  p  ‘-—A  X  =  4p  ‘  a(pp  ‘ )  x  =  -j^(p  'ap  )(p  '  X  ) 


We  want  to  get  Equation  (2.1 )  into  a  simpler  form  to  make  it  easier 


-1 


to  solve.  To  do  this  we  will  let  U  =  P  X  .  This  matrix  equation  can 
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-1 

be  easily  differentiated  with  respect  to  time,  since  P  is  a  constant 


matrix.  The  first  derivative  with  respect  to  time  is  U  =  P  X  . 

Since  Equation  (2. 1 )  is  a  second  order  differential  equation,  taking  a 

second  derivative  with  respect  to  time  yields  U  =  P  X  .  We 

introduce  the  vector  variable  U  into  Equation  (2.1)  by  substituting 
-¥ 

U  =  P”*  X  into  the  left  side  of  this  equation.  Then,  if  we  substitute 

- 1  ^  _i 

P  AP  =  D  and  U  =  P  X  into  the  right  side.  Equation  (2.1)  becomes 
^  k 

U  =  -j^D  U  .  This  system  of  second  order  differential  equations  is 


easier  to  solve  than  X  =  A  X  . 


Exercise  2.1 


Show  Why  the  system  U  =  '  0  _,]  =-^D 


U  is 


easier  to  solve  than  X 


X.  1  -2  J 


•  — X 

"  m  ^  ^ 


Recall  that  we  are  trying  to  solve  for  the  vector  X  . 

Rearranging  U  =  P  X  ,  we  get  X  =  P  U .  which  tells  us  that  instead 
— ^  — > 

of  finding  X  we  need  only  find  PU.  Since  we  already  know  the 
matrix  P,  we  must  find  the  vector  U  .  If  we  multiply  both  sides  of 


the  matrix  equation  U  =  — D  U  by  the  matrix  P,  we  get 
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P  U  =  -^P  D  U .  Rewriting  the  left  side  of  p  U  =  -^P  D  U .  using  the 


m 


.(i) 


notation  where  P  represents  the  i  th  column  (i  -1  and  2)  of  the 
matrix  P,  we  obtain 


P  U 


(p(l)  p(2)) 

.“2. 

-  n<2)  r. 

=  P  U  P 


Rewriting  the  right  side  of  P  U  =  -^P  D  U .  we  obtain 

^Pd3=^[p‘'\u,h-p'%uJ 


Thus  P  U  =^P  D  U  can  be  written  as 


ji)  -  „(2) .  i 


We  can  simplify  this  equation  by  multiplying  through  by 


m 


gathering  terms  and  moving  all  terms  to  the  left  side. 

[pf  ■>  a ,  -  i  p"’>.,  u ,]  ^p'^' 0  2  p'%  u  2]  =  0 


r- 

Factoring  out  P  and  P  .  we  have 
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Since  the  columns  of  P  are  eigenvectors  of  A  which  correspond  to 
distinct  eigenvalues,  we  know  they  are  linearly  independent  (in  fact 
they  are  orthogonal).  The  equation  produces  a  finite  linear 
combination  of  linearly  independent  vectors  which  equals  zero,  thus 
the  coefficients  of  ^  and  P  must  be  zero.  If  we  set  each  of  the 


coefficients  in  the  equation  above  equal  to  zero,  we  obtain 


These  are  both  second  order  linear  differential  equations  which  can 
be  solved  using  basic  techniques.  (See  Appendix  A.)  If  we  let 

r.  =*-  where  i=  1 ,  2.  then  these  equations  become 


u,  +  r.u.  =  0  and  0.,+ r.,  u  =  0 
III  222 


Using  the  following  formulas,  we  can  solve  for  the  vector  U  .  (Note: 
To  determine  whether  rj  is  zero,  negative  or  positive,  substitute 

into  r.  =- -^X  .) 

I  ni  1 
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If  r.  =  0. 
1 

If  r.<0, 

I 

If  r.  >  0, 


then  u.  =  c., 

I  il 


then  u.  =c.. 

i  il 


then  u.  =  c., 

I  il 


t  +c 


i2 


e 


cos 


+  c..e 

i2 


sin 


The  solution  to  the  original  system  of  differential  equations 
■ft  k 

X  =  A  X  is  found  by  substituting  the  values  for  both  the  matrix  P 

— >  — > 

and  the  vector  U  into  the  equation  X  =  P  U  . 


Exercise  2.2 

Using  the  above  technique,  solve  the  following  system  of  differential 


equations  X  = A  X  .  Where  X  = 


X  = 


^  X 


.  That  is.  find  the  two  equations  which  describe  X|  and  X2. 


Exercise  2.3 

Given  a  horizontal  spring-mass  system  similar  to  Figure  2.3  with  n 
blocks,  the  following  equations  would  express  the  kinetic  and 
potential  energies  of  the  system. 


Kinetic  Energy 


Potential  Energy 
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I  ”  7 

T  =  y  m  X. 

^  i  =  1  ' 

+  2b,2X,X2  +  2b,3X,X3+...+  2b,_,x,x_, 

+  2b33X3X3  +  2b3^X3X^+...+  2b3_^X3X„ 


+  2b.  .  .  x.x.  .  + ...  +  2b.  X.  X 
I.  i+l  I  1+1  i.n  1  n 


+  2b  .X  ^x  .+2b  ^  X  .,x 

n-2.  n-I  n-2  n-l  n-2,  n  n-2  n 


+  2b  ,  X  X  1 
n-l,  n  n-l  n 


Use  the  equations  of  motion  to  find  the  system  of  differential 
equations  which  model  the  spring-mass  system  with  n  blocks.  To 
solve  this  system,  generalize  the  procedure  used  to  solve  the  system 
of  differential  equations  which  model  the  spring-mass  system  with 
two  blocks.  (Hint:  Some  of  the  material  that  has  been  discussed  can 
be  used  directly,  while  other  portions  will  need  some  modifications.) 


Another  aspect  of  the  spring-mass  system  that  we  want  to 
consider  is  the  oscillations  of  the  system  as  a  whole.  From  our  work 
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above,  we  know  the  solution  to  the  second  order  differential 

Tt  j.  -»  ->  -» 

equation  X  =  A  X  can  be  found  by  using  X  =  P  U  .  In  Exercise 
2.2  we  found 


X|  =Uj+U2  and  X2  =  -Uj+U2. 

This  is  a  system  of  two  linear  equations  which  we  can  solve  for  u  j 
and  U2-  Thus  we  have  the  equations 


(2.2) 


and  0^  = 


each  of  which  gives  a  relationship  between  the  variables  Xj  and  X2. 
It  is  important  to  note  that  we  could  have  found  these  equations 

_i 

directly  from  the  matrix  equation  U  =  P  X  ,  but  this  would  involve 
finding  P  V  Using  Figure  2.5,  we  can  recall  the  configuration  of 

this  spring-mass  system.  Since  the  springs  were  stretched  by  • 
differing  amounts,  a  different  frequency  (the  number  of  vibrations 
per  unit  time)  is  associated  with  each  of  the  variables  Xj  and  X2. 
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-A/W - 1 


I - Xj - 1 

Frequency  1 


hX2H 

Frequency  2 


Figure  2.5 


This  spring-mass  system  has  two  separate  modes  in  which  it 
vibrates. 

X  -X2  x,-X2 

In  the  first  mode  =  — 2 —  U2=0.  Since  — 2 — 

represents  the  how  the  distance  between  the  two  blocks  is 
changing,  the  first  mode  of  vibration  describes  how  the  distance 
between  the  two  blocks  is  changing.  For  instance  if  X2  is  greater 

than  Xj,  then  the  change  in  distance  between  the  blocks  is  smaller 

than  the  distance  between  the  blocks  when  the  spring-mass  system 
is  in  equilibrium.  However,  if  Xj  is  greater  than  X2.  then  the  change 

in  distance  between  the  blocks  is  larger  than  the  distance  between 
the  blocks  when  the  spring-mass  system  is  in  equilibrium.  Thus, 
the  oscillation  of  the  system  in  this  mode  is  described  by  how  the 
distance  between  the  two  blocks  is  changing  which  corresponds  to 
the  frequency  associated  with  the  second  eigenvalue  ^2.  To 

visualize  this,  consider  the  series  of  “snapshots"  of  the  spring-mass 
system  in  motion  in  Figure  2.6,  where  the  banner  is  made  of  an 
elastic  material  and  indicates  the  distance  between  the  two  blocks 
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in  the  system.  When  this  system  vibrates,  we  would  see  the  banner 


/V/WW- 


^ 1 

Figure  2.6 


^  2  12 

In  the  second  mode  u  j  =0  and  0^=  — 2 —  Since  — 

represents  how  the  center  of  mass  of  the  system  has  changed,  the 
second  mode  of  vibrations  describes  the  displacement  of  the  center 
of  gravity.  Thus  the  oscillation  of  the  system  in  this  mode  is  where 
the  center  of  mass  of  the  system  vibrates  at  the  frequency 
associated  with  the  first  eigenvalue  Xj.  To  visualize  this,  consider 

the  series  of  diagrams  in  Figure  2.7.  where  the  flag  indicates  the 
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center  of  mass  of  the  system.  When  this  system  vibrates,  we  see  the 
flag  moving  back  and  forth  with  a  frequency  associated  with  Xj. 

This  is  indicated  by  the  following  series  of  “snapshots”  of  the  spring- 
mass  system  in  motion. 


Exercise  2.4 

Suppose  the  spring-mass  system  we  have  been  studying  was  lying 
free  in  the  xy-plane,  that  is,  the  ends  of  the  springs  are  not 
anchored.  Figure  2.8  can  help  us  visualize  this. 


7-axis 
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./W/S/W- 


x-axis 


Figure  2.8 


Using  the  information  we  have  gained  by  studying  the  stationary 
spring-mass  system,  describe  the  motion  (including  the  vibrations) 
that  can  occur.  Note,  there  is  no  need  to  find  the  frequencies  to 
complete  this  exercise.  (Hint:  consider  other  types  of  motion, 
besides  vibrations.) 


Exercise  2.5 

(a)  Determine  the  system  of  differential  equations  that  model  the 
motions  of  the  spring-mass  system  given  in  Figure  2.9. 

(b)  Solve  the  system  of  differential  equations. 

(c)  Describe  the  possible  configurations  in  which  it  vibrates. 


Sping  1  m  Spring  2  m  Spring  3  m  Spring  4 


Figure  2.9 
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Exercise  2.6 

Suppose  the  spring-mass  system  of  Exercise  2.5  was  lying  free  in 
the  xy-plane,  that  is,  the  ends  of  the  springs  are  not  anchored. 
Describe  the  motion  (including  vibrations)  that  can  occur  for  this 
system.  Compare  these  motions  with  the  motions  found  in  Exercise 
2.4. 


Section  3  A  Closed  Spring-Mass  System 

In  this  section  we  will  discuss  how  to  mathematically  model 
the  spring-mass  system  in  Figure  3.1  and  determine  the  possible 
motions  of  the  system.  This  system  lies  in  the  xy-plane  with  none  of 
its  blocks  anchored.  The  mass  of  each  of  the  three  blocks  is  the 
same  and  is  denoted  by  m.  L  is  the  length  of  each  spring  when  the 
system  is  in  its  equilibrium  configuration  and  k  is  the  spring 
constant,  which  is  the  same  for  each  spring. 
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^2 


Figure  3.1 

This  system  is  in  stable  equilibrium  when  the  displacements 

^1  “y  energy  of  the  system  in  Figure 

3.1  we  need  to  find  the  kinetic  and  potential  energies  of  the  system. 
Recall  that  the  kinetic  energy  of  the  system  is  one  half  the  mass 
times  the  sum  of  the  square  of  the  first  derivative  of  each  of  the  six 
variables  with  respect  to  time.  Thus,  the  kinetic  energy  is 

_  I  f  .2  .2  .2  .2  .2  .2 

T=2m(x,  +  X24X3  +  y,  +  y2  +  y3 

Finding  the  potential  energy  requires  more  work.  Since  the 
potential  energy  of  the  system  is  the  sum  of  the  potential  energies  of 
the  springs,  we  first  need  to  find  the  potential  energy  of  each  spring. 
We  will  consider  each  side  of  the  triangle  individually. 

The  first  side  that  we  look  at  is  given  in  Figure  3.2. 
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Figure  3.2 

The  potential  energy  of  this  spring  is  one  half  the  spring  constant  k 
times  the  square  of  the  distance  that  the  spring  is  stretched.  If  we 
let  d  represent  the  length  of  the  spring  after  it  has  been  stretched. 

then  the  displacement  of  the  spring  from  its  equilibrium  position 
(the  distance  that  the  spring  is  stretched)  is  |d-L|.  Expressing  the 

potential  energy  for  the  spring  in  terms  of  |d-L|.  we  have 

V,2  =  lk|d-L|' 

Where  the  subscript  12  of  V  indicates  that  we  are  finding  the 
potential  energy  of  the  spring  that  is  stretched  between  the  block 
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with  coordinates  Xj  and  yj  to  the  block  with  coordinates  X2  and  y2. 
Now.  we  want  to  rewrite  V  j  2  using  the  variables  x  j .  y  j ,  X2.  and  y2- 
To  do  this,  we  must  first  simplify  the  expression  for  the  distance  d. 

1 

1 

=  ■.[2  ) 

r  2  2 

=  |tL'  +  l(x,-xJ  +  (x,-xJ  H.iL'-V3L(y,-yJ+(y,-yJ  • 

r  ^ 
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using  its  Taylor  series  expansion.  The  terms  of  higher  powers  have 
been  grouped  together  for  convenience. 


d=  L 


1  + 


2  [("  I  - 


1  ! 


+  terms  of  higher  powers 


Exercise  3.1 

Verify  that  the  expression  above  is  indeed  the  Taylor  series 

expansion  for  the  quantity  (shown  below).  Hint:  write  f 

in  its  Taylor  series  expanded  about  zero.  Recall  that  L  is  the  length 
of  the  spring  in  equilibrium,  thus  L*0.  (Hint:  To  make  it  easier  to 

take  the  derivative  of  f  let  r=  ^  derivative  of  f(r).) 


The  expression  preceding  Exercise  3.1  can  be  simplified  by 
multiplying  through  by  L  and  then  moving  L  to  the  left  side.  The 
resulting  quantity  is  what  we  want. 
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d-  L  =  -j[^(^X|-  ^2)“  V3(y  >^2)]+  terms  of  higher  powers 

This  quantity  can  now  be  substituted  into  the  formula  for  the 
potential  energy  V  j  2- 

V,2  =  7k|d-L|' 

=  -x^]- VT^y^ -y^^  +  terms  of  higher powersj  j 

I  =Y|-^[^Xj-X2]-V3(^yj-yJj^+ terms  of  higher  powers  | 

Recall  from  an  earlier  discussion  that  the  Taylor  series  expansion  for 
potential  energy  can  not  have  any  nonzero  linear  terms  because  we 
are  in  a  system  which  has  an  equilibrium  configuration.  Also,  we 
are  only  considering  small  vibrations  so  we  ignore  the  terms  of 
higher  powers.  The  formula  for  potential  energy  V  j  2  is 


Jc 

2 


22  22 

^  ''2  1 

1.4  ^  4  ^  — ^  —  -2"l"2-— "1^1^— "2^1 


VI . VI 


VI  VI 

'*'~2  ^  V2 - ^^2^2 


3 


Exercise  3.2 

Using  Figure  3.3,  find  V23.  Hint:  the  procedure  is  similar  to  the  one 
used  to  find  V  j  2 


Figure  3.3 

Exercise  3.3 

Using  Figure  3.4,  find  Vj^. 


Figure  3.4 


As  stated  before,  the  potential  energy  of  the  system  is  the  sum  of  the 
potential  energy  of  each  spring.  Thus  we  have 
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Since  this  system  requires  six  coordinates  to  fully  describe  it,  we 
know  we  must  have  six  equations  of  motion.  These  are 


First,  we  find  the  left  side  of  each  equation  of  motion,  then  find  the 
right  side  and  equate  the  two.  Thus,  we  have  the  following  six 
equations. 
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=  m  { 


5  1  VI  VI  1 

^x,  +  -X2  +  X34.— y,-— y2| 


Y  -  JL  I  ±Y  _  1y  4-  -Lv  _  ^  I 

^2"m  l4  I  2  2  4  3  4  4  ^3j 


k  I  1  5 

^3“  m  4^2  4 


X  *:^y  -:^y 

^3  4  ^2  4  ^3 


} 


y,=  rnl  — 


k  /VI.,  VI  _  3.  3.  1 

*2  4^1  4^2/ 


v>  k  /  VT..  ,  VI  ■■  .  _3,,  3.,,  .3.^  I 
^2  ml  4  4  ^3"^  4^1  2^2'^4^3j 


V 

m  1  4  2  4  3 


y,=  7^1-^X,,-^X,+  ly^-lyj 


This  system  of  six  equations  can  be  written  as  a  matrix  equation.  In 

order  to  eliminate  fractions  from  the  matrix,  we  factor  ttt  out  of 

4m 

each  equation,  which  results  in  being  factored  out  of  the 
coefficient  matrix  A. 


^  -5 

1 

4 

VI 

-VI 

0  ' 

1 

^2 

1 

-2 

1 

-VI 

0 

VI 

^2 

"3 

k 

4 

1 

-5 

0 

VI 

VI 

y, 

4m 

VI 

-VI 

0 

-3 

3 

0 

yi 

y^ 

-VI 

0 

VI 

3 

-6 

3 

y^ 

2 

U3J 

k  0 

VI 

-VI 

0 

3 

-3  > 

2 

U3J 
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The  generalized  theory  developed  in  Exercise  2.3  describes  the 
situation  when  n=6.  Thus,  we  begin  by  finding  the  eigenvalues  and 
eigenvectors  of  the  symmetric  matrix  A.  One  way  we  could  proceed 
would  be  to  use  the  sixth  degree  characteristic  polynomial  to  find 
the  eigenvalues  directly.  However,  this  would  require  finding  the 
determinant  of  a  6x6  matrix.  Using  the  cofactor  expansion  method 
would  require  6!  or  720  calculations  to  find  the  value  of  the 
determinant.  We  could  also  use  a  computer  program.  For  example, 
the  user's  guide  to  the  computer  program  LINPACK  (Dongarra, 
Bunch  and  Stewart,  1979)  describes  how  the  program  can  be  used 
to  approximate  the  eigenvalues  and  eigenvectors  of  the 
characteristic  polynomial.  This  would  be  quicker,  but  would  not 
give  us  any  insight  into  the  possible  types  of  vibrations  of  the 
system.  Instead,  let  us  consider  the  symmetric  matrix  A  and  see  if 
we  can  use  our  knowledge  of  matrices  to  reduce  the  amount  of 
work  required  to  find  the  eigenvalues.  In  general,  the  coefficient 
matrix  which  represents  an  application  is  much  larger  than  a  6  x  6 
matrix,  but  is  still  a  symmetric  matrix.  The  approach  used  by 
applied  mathematicians  working  on  large  systems  would  be  tO:  1) 
manually  work  through  the  theory  of  a  smaller,  related,  and  less 
complicated  system,  2)  enlarge  the  system  and  use  a  computer  to 
find  the  eigenvalues  and  eigenvectors,  3)  interpret  the  physical 
meaning  of  the  information  from  the  computer  by  comparing  the 
results  with  the  results  found  in  step  1,  and  finally,  4)  change  the 
model  so  that  it  reflects  the  desired  system  as  closely  as  possible. 

For  example,  in  a  more  complicated  system  not  all  of  the  blocks  may 
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be  of  the  same  mass,  nor  the  springs  be  of  the  same  length  or  have 
the  same  spring  constant.  Step  1  may  be  to  consider  a  system 
where  all  of  the  blocks  have  the  same  mass,  the  springs  are  all  of  the 
same  length,  and  each  spring  has  the  same  spring  constant. 

Therefore,  we  will  start  our  work  by  finding  the  determinant 
of  the  matrix  A-XI 


-5-1 

1 

4 

VI 

-VI 

0 

1 

-2-1 

1 

-VI 

0 

VI 

det(A-  11)  = 

4 

1 

-5-1 

0 

VI 

-a/I 

VI 

-VT 

0 

-3-1 

3 

0 

-VI 

0 

VI 

3 

-6-1 

3 

0 

VI 

-VI 

0 

3 

-3-1 

Next,  we  replace  the  first  row  by  the  sum  the  first  three  rows  and 

replace  the  last  row  by  the  sum  of  the  last  three  rows  to  obtain  the 

following  interesting  matrix. 

-X 

-1 

-1 

0 

0 

0 

1 

-2-1 

1 

-VI 

0 

VI 

4 

1 

-5-1 

0 

VI 

-VI 

I  a-  I 

VI 

-VI 

0 

-3-1 

3 

0 

=  iA  1 

-VI 

0 

VJ 

3 

-6-1 

3 

0 

0 

0 

-1 

-1 

-1 

If  X  were  set  equal  to  zero,  the  matrix  A',  defined  above,  would  have 
two  rows  of  zeros  indicating  that  X=0  is  an  eigenvalue  of  A  with 
multiplicity  ^  2. 
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We  pause  for  a  moment  in  our  pursuit  of  eigenvalues  to  find 
the  eigenvectors  associated  with  X=0.  To  begin,  A,=0  is  substituted 
into  A'  so  that  the  matrix  equation  A‘  X  =  0  ,  which  has  been  written 
in  augmented  form,  can  be  solved. 


"  0 

0 

0 

0 

0 

0 

1 

0 " 

1 

-2 

1 

-VI 
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VI 
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0 
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1 

-5 
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-VI 
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-VI 
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-3 

3 

0 

1 

0 

-VJ 

0 

VI 

3 

-6 

3 

1 

0 

.  0 

0 

0 

0 

0 

0 

1 

0  > 

Using  Gaussian  elimination,  we  reduce  this  system  to  a  form  that 
can  easily  be  solved. 


^  1 

0 

-  1 
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0 

.  0 

0 

0 

0 

0 

0 

0  > 

From  this  augmented  system  of  equations,  with  three  rows  of  zeros, 
we  know  three  of  the  variables  (x^. 

value,  forcing  the  remaining  three  variable  (x  j ,  X2  and  y  j )  to  take 


y2  and  y^)  can  take  on  any 


55 


on  specific  values  given  by  the  following  equations,  obtained  from 
the  augmented  matrix  in  Equation  (3.1).. 

(3.2)  X2=X3-V3y2  +  V3y3 

y,  =  2y2-y3 

Thus,  by  letting  x^,  y2  ^nd  y^  take  on  specific  values,  we  will  have 

three  linearly  independent  eigenvectors.  This  means  the  eigenvalue 
X=0  must  have  multiplicity  three. 

Before  we  actually  determine  the  values  of  the  eigenvectors, 
let  us  pause  for  a  moment  to  see  how  we  can  rewrite  the  potential 
energy  function  in  a  slightly  different  format  which  will  help  us  to 
determine  its  value  under  certain  conditions.  Recall,  the  potential 
energy  of  the  system  is  the  sum  of  the  potential  energy  of  each 
spring.  Thus  we  have 
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We  are  able  to  rewrite  the  potential  energy  in  this  format  because  A 
is  a  symmetric  matrix.  We  define  a  function  which  can  be  rewritten 
in  this  fashion  as  a  quadratic  form.  Thus,  Equation  (3.3)  is  the 
potential  energy  expressed  as  a  quadratic  form. 


(3.3)  v(D  =  X^^AX 

8m 

By  the  definition  of  an  eigenvector  X  ,  which  is  associated 
with  the  eigenvalue  X  of  the  matrix  A.  we  know  that  A  X  =X  X  .  If 
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we  let  the  three  eigenvectors,  associated  with  the  eigenvalue  X.=0,  be 


represented  by  X,  „  •  X  .and  X  then 

^  X  =0  X  =0  X  =0 

1  2  3 


T  T  T 

V  X,  ,  =x,  ,t^AX,  ^=X,  ,=  x, 

X  =0  X  =0  8  m  X  =0  X  =0  8  m  x  =o  x  =o  8m  x 

{  j  y  J  J  j  j  J  j 


for  j=  1 .  2  or  3.  This  tells  us  that  the  potential  energy  is  zero.  We 
now  examine  the  physical  interpretation  of  zero  potential  energy. 

To  have  zero  potential  energy  in  the  system,  all  the  springs 
must  remain  the  same  length  L  as  in  equilibrium.  Thus,  the  only 
type  of  motion  possible  occurs  when  the  entire  system  moves  as  a 
unit.  This  is  called  a  rigid  motion.  Since  the  spring-mass  system 
lies  in  the  xy-plane,  there  are  only  two  types  of  rigid  motion: 
translations  (movement  in  the  x-  or  y-direction  only)  and  rotations 
(the  system  pivots  around  its  center  of  mass).  These  two  motions 
can  also  be  combined. 

— > 

If  we  consider  the  vector  X  .  as  a  translation  in  the  x-direction 
only,  then  the  variables  x  j .  X2  and  x^  must  all  change  by  the  same 

value  and  the  variables  y  1 .  y2  not  change.  We  can 
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express  this  vector  as 


c 

c 

0 

0 

.  0  > 


If  we  let  c=  1 .  an  eigenvector  associated 


with  the  eigenvalue  X=0  is  X 

X  =0 
1 


^  1 
1 
1 

0 

0 

.0 


.  The  graphical 


J 


interpretation  of  ^  q  can  be  seen  in 
Figure  3.5. 


Xj-l 


Resulting 

vector 


Figure  3.5 
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If  the  translation  is  to  the  right  (in  the  positive  x-direction).  then  oO 
and  if  it  is  to  the  left  (in  the  negative  x-direction),  then  c<0. 


Exercise  3.4 


Verify  that  X 


X  =0 

1 


^  1  ^ 

1 

1 

0 

0 


is  a  solution  to 


VOJ 


X  =  X  '  ^ 

X  =0  X  =0  8  m 
J  /  J 


AX  =  0 
X  =0 
j 


The  other  type  of  translational  motion  we  wish  to  consider 
occurs  when  the  system  moves  in  the  y-direction  only.  The 
variables  x  j ,  X2  and  x^  do  not  change  while  the  variables  y  j .  y2  2nd 

y^  must  all  change  by  the  same  value.  We  can  express  this  vector 


(0] 

0 

0 

c 

c 

.  C  > 


as 


If  we  let  c=  1 ,  an  eigenvector  associated  with  the 


60 


eigenvalue  X-0  is 


x^=o 


f  ^  1 

0 

0 

I 

1 

.  I  > 


The  graphical  interpretation  of 


X  ,  can  be  seen  in 
Figure  3.6. 


System  Motion: 

Translation  along 


Resulting 

vector 


If  the  translation  's  upward  (in  the  positive  y-direction),  then  oO 
and  if  the  translation  is  downward  (in  the  negative  y-direction), 
then  c<0. 
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Exercise  3.5 


Verify  that  X 


X  =0 
2 


0 
0 
1 
1 

1 ; 


is  a  solution  to 


X  =0 
V  j  y 


~*T  —  k 

=  X,  A  X  =0. 

^  =0  8m  X  =0 
J  j 


Taking  a  careful  look  at  the  vectors  X,  and  X  we  see 
®  X  =0  X  =0 

I  2 

they  are  orthogonal.  That  is.  their  dot  product  is  zero.  Since  the 
eigenspace  associated  with  the  eigenvalue  X=0  has  dimension  at 
least  three,  we  know  there  is  a  third  linearly  independent 
eigenvector  associated  with  the  eigenvalue  X=0.  There  are  two 
ways  we  could  proceed  at  this  point.  The  first  is  to  use  the  three 
equations  in  (3.2)  and  choose  values  for  x^.  y2  and  y^-  For 

example,  let  x^'O,  y2=l  and  y3=0.  then  use  the  Gram-Schmidt 
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process  to  find  a  vector  which  is  orthogonal  to  both  ^  and 

A  "U 

I 


X  .  The  other  way  is  to  replace  two  of  the  rows  of  zeros  m  the 
X  =0 
2 


coefficient  matrix  in  Equation  (3.1 )  by  the  eigenvectors  X  and 

X  —0 
\ 


X  .  The  solution  to  this  new  augmented  matrix  must  satisfy  all 
2~ 

the  equations  which  form  the  augmented  matrix.  Hence,  the 
solution  to  the  augmented  system  will  satisfy  both  Xj  ♦X2*X3=0  (from 


X,  )  and  y, +y-,*yT=0  (from  X  ).  A  vector  whose  entries 
A.  =0  1  2  3  X,  =0 

i  2 

satisfy  both  of  these  equations  is  orthogonal  to  the  eigenvectors 

X  and  X  .  Also,  from  these  two  equations,  we  see  in  the 
X  =0  X  =0 

I  2 

solution  to  the  equations  associated  with  the  augmented  matrix,  the 
Xj  values  must  sum  to  zero.  Therefore,  there  is  no  translational 

motion  in  the  x-direction.  Similarly,  there  is  no  translational  motion 
in  the  y-direction.  Thus,  the  center  of  mass  does  not  move.  Since 
this  motion  is  a  rigid  motion  (X=0)  and  the  center  of  mass  of  the 
system  does  not  move,  the  rigid  motion  must  be  a  rotation. 


Reducing  the  following  augmented  matrix  which  is  EquationO.I ) 


T  T 

->  -» 

with  two  of  its  rows  of  zeros  replaced  by  X  and  X 


^  1 

0 

- 1 

0 

0 

0 

1 

0  ^ 

0 

1 

- 1 

0 

VT 

-VI 

1 

0 

0 

0 

0 

I 

-2 

1 

1 

0 

1 

1 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

1  0 

0 

0 

0 

0 

0 

1 

0  > 

we  obtain 


f 

1 

0 

0 

0 

0 

VJ 

3 

\ 

0 

0 

1 

0 

0 

0 

2V3 

3 

0 

0 

0 

I 

0 

0 

VT 

3 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

<  0 

0 

0 

0 

0 

0 

0  > 

This  augmented  matrix  above,  can  be  interpreted  as  the  following 
equations. 


X 


3 


X 


VI 

— y 


3 


y 


y 


If  we  let  y3=3,  then  the  resulting  eigenvector  is 


S=o 


'  -VI ' 
2V3 
-VI 

-3 

0 

<  3 


The  graphical  interpretation  of 


can  be  seen  in  Figure  3.7. 
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Exercise  3.6 


Verify  that  X 


X  =0 
3 


^  -  VI  ^ 
2V3 
-VI 

-3 
0 
3 


is  a  solution  to 


V 


X  =0 

V  j  J 


=  X 


T  -k 


AX  =  0 


^  =0  8  m  A,  =0 
J  j 
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Also,  verify  that  X  is  orthogonal  to  both  X ,  and  X  ^ 

^  X.  =0  X  =0  X  =0 


So  far  we  have  found  only  three  eigenvalues  and  their 
associated  eigenvectors.  The  remaining  three  eigenvalues  can  be 
found  using  the  determinant  of  the  matrix  A-XI  which  can  be 
reduced  to  |a  |.  For  convenience.  jA’]  has  been  repeated  below. 


-X 

-X 

-X 

0 

0 

0 

1 

-2-X 

1 

-VI 

0 

VI 

4 

1 

-5-X 

0 

VI 

-vT 

VI 

-VI 

0 

-3-X 

3 

0 

-V3 

0 

VI 

3 

-6-X 

3 

0 

Q 

0 

-X 

-X 

-X 

Using  Gaussian  elimination,  we  will  reduce  the  matrix  to  a  form 
which  will  make  the  determinant  easier  to  find.  We  will  use  only 
row  (or  column)  operations  that  do  not  change  the  value  of  the 


determinant. 

After  several  row  operations. 

-  — 

we  obtain 

1 

1 

1 

0 

0 

0 

0 

-3-X 

0 

-2VI 

-VI 

0 

0 

-3 

-9-X 

VI 

2VI 

0 

0 

-2VI 

-VI 

-3-X 

3 

0 

0 

-3VI 

0 

-6-2X 

-3-X 

0 

0. 

0 

0 

1 

1 

1 
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At  this  point,  we  could  find  the  determinant  using  the  cofactor 
expansion  method.  However,  if  we  do  one  column  operation  we  will 
greatly  reduce  the  number  of  calculations  needed.  We  add  -2  times 
the  fifth  column  to  the  fourth  column  producing  a  new  fourth 
column. 


1 

1 

1 

0 

0 

0 

0 

-3-X 

0 

0 

-VI 

0 

0 

-3 

-9-X 

-3  VI 

2VI 

0 

0 

-2VI 

-VI 

-9-X 

3 

0 

0 

-3VI 

0 

0 

-3-X 

0 

0 

0 

0 

1 

1 

1 

We  are  now  ready  to  use  the  cofactor  expansion  method  to  find  the 
determinant  of  the  matrix  A'.  Expanding  by  the  first  column  we 
have 


A'  I  =  1  (-  1) 


1+1 


-3-X 
-  3 

-2V3 

-3V3 

0 


0 

-9-X 

-VI 

0 

0 


0  -VI  0 
-  3VI  2VI  0 

-9-X  3  0 

0  -3-X  0 

1  11 


Now,  expand  the  resulting  cofactor  by  the  fifth  column.and  obtain 


-3-X 

0 

0 

-VJ 

I-  K-i)^^^ 

-3 

-2VI 

-9-X 

-VI 

-3VI 

-9-X 

2V3 

3 

-3VI 

0 

0 

-3-X 
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Next,  expand  the  resulting  cofactor  by  the  first  row  .and  obtain 


A'  I  =  1  •  1 


(_3-X)(-  I) 


1+1 


+  -  V3(-  1) 


1+4 


9-X  - 

-3V3 

2VI 

-9-X 

3 

0 

0 

-3-X 

-3 

-9-X 

-3  VI 

-2V3 

-VI 

-9-X 

-3  VI 

0 

.  0 

=  [(-9-X)  ^-9][(-3-X)^-9]. 


To  find  the  eigenvalues  of  the  original  matrix,  we  set  each  factor 
equal  to  zero  and  solve  for  X. 

(-9-X)^-9  =  0  (-3-X)^-9=0 

(_g_X)^=9  (-3-X)^  =  9 

-9-X=±3  -3-X=±3 

X  =  -  12,-6  X  =  -6,-  0 

Since  we  have  already  determined  that  the  eigenvalue  X=0  has 
multiplicity  of  at  least  three,  the  fact  X=0  occurs  above  should  be  no 
surprise.  The  remaining  eigenvalues  for  the  matrix  A  are  X=-12  and 
X=-6,  the  latter  with  multiplicity  2. 
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To  help  us  find  the  associated  eigenvectors  for  the  remaining 

— ►  “4 

three  eigenvalues,  we  recall  the  equation  A  X  =A.  X  from  the 

Ar  A 

definition  of  an  eigenvector  X  ,  and  Equation  3.3 

V  X  =X  — A  X  .  which  describes  the  potential  energy  as  a 
\  X  y  X  8  m  X 

quadratic  form  using  eigenvectors.  As  we  saw  earlier,  these  two 

equations  can  be  combined  as  V  X  =  X  — ^  X  X  .  We 

\  X  y  X  8  m  X 

observe  that  the  only  way  this  equation  can  equal  zero  is  if  X=0  or 


X  ^  is  the  zero  vector.  However,  since  we  are  only  looking  at  X=- 1 2 


or  X»-6,  which  are  nonzero  values,  we  must  have  that  X  be  the 

X 

zero  vector  in  order  for  V^X  ^  ^  to  equal  zero.  Clearly  this  cannot 


happen  because  X  is  an  eigenvector  which  by  definition  is  never 

X 

equal  to  the  zero  vector.  This  indicates  that  the  potential  energy  of 
the  system  is  not  zero.  Hence,  the  potential  energy  of  each  spring  is 
not  zero,  so  the  length  of  at  least  one  of  the  springs  must  change. 
Thus,  we  do  not  have  a  rigid  motion.  Also,  we  recall  that  the 
determinant  of  the  matrix  A-XI  can  be  reduced  by  summing  the  first 
three  rows  and  the  last  three  rows  to  give  I  A’l. 
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-X 

-X 

-X 

0 

0 

0 

1 

-2-X 

1 

-V3 

0 

VI 

4 

1 

-5-X 

0 

V3 

-VI 

V3 

-V3 

0 

-3-X 

3 

0 

-VJ 

0 

VJ 

3 

-6-X 

3 

0 

0 

0 

-X 

-X 

-X 

If  we  substitute  in  X=-12  or  X=-6,  the  first  row  will  contain  constant 
values  for  Xj.  X2.  and  x^  and  the  last  row  will  contain  constant 

values  for  Vj,  ^2'  >^3  Ffoni  an  earlier  discussion  (following 

Exercise  3.5)  this  indicates  there  is  no  translational  motion  in  either 

the  X-  or  y-directions,  so  we  know  the  center  of  mass  does  not  move. 

Thus  the  motion  associated  with  the  last  three  eigenvectors  can  be 

thought  of  as  vibrations  of  the  blocks  (but  not  a  translation  or 

rotation)  with  the  center  of  mass  remaining  fixed. 

First,  we  find  the  two  linearly  independent  eigenvectors 

associated  with  the  eigenvalue  X=-6.  Since  X=-6  has  multiplicity 

— > 

two,  the  solution  space  of  the  augmented  matrix,  (A-XI)  X  =  0  or 
(A +  6  0  X  =  0  will  have  dimension  four.  That  is,  when  the 
augmented  system  is  reduced,  we  will  have  two  rows  of  zeros. 

Thus,  four  of  the  variables  can  be  written  in  terms  of  two  of  the 
other  variables.  These  two  variables  can  be  assigned  values  which 
will  produce  two  linearly  independent  eigenvectors.  If  we  let  X2 


and  y2  be  these  two  variables,  then  Xj,  yj.  x^  and  y^  can  be  written 
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in  terms  of  X2  and  ^2-  assign  values  to  X2  and  ^2 

assured  of  getting  a  linearly  independent  eigenvector,  is  to  first  let 
X2=0  and  y2=  1 .  and  then  let  X2=  I  and  y2=0-  Let  us  consider  the 

geometric  interpretation  of  these  cases. 

CASE  1.  X2=0andy2=l 

Since  the  center  of  mass  for  this  configuration  remains  fixed,  the  ^2 
component  must  be  balanced  by  the  sum  of  the  yj  and  the  y-^ 

components.  Because  X2=0.  we  know  the  components  Xj  and  x^ 

must  be  equal  in  magnitude  and  of  opposite  sign.  These  components 
can  be  seen  in  Figure  3.8. 
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System  Motion: 

Stationary  Vibration 


Resulting 

vector 

72-1 


Resulting 

vector 


Figure  3.8 


Letting  c=l,  one  eigenvector  associated  with  X=-6  is 


^  1  ^ 
0 

-  1 


CASE  2.  X2=I  and  y2  =  0 
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Since  the  center  of  mass  does  not  move,  the  X2  component  is 
balanced  by  the  sum  of  the  Xj  and  x^  components.  Because  y2=0. 

we  know  the  components  yj  and  y^  must  be  equal  in  magnitude 

and  of  opposite  sign.  These  components  can  be  seen  in  Figure  3.9. 


Figure  3.9 

Letting  c=l,  a  second  eigenvector  associated  with  A,=-6  is 

1  ^ 

2 
1 

j_ 

2  ■ 

1 

0 

1  > 


System  Motion:  ^2  “  ® 
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Exercise  3.7 

— > 

Show  that  X  and  X  are  orthogonal  (their  dot  product  is 
X|=-6 

zero). 


Therefore.  {X  ,  X  }  is  a  set  of  orthogonal  eigenvectors 

A.  ^=-0  A,  ^=-6 

associated  with  the  eigenvalue  X,=-6. 

It  remains  for  us  to  find  the  single  eigenvector  associated  with 
the  eigenvalue  X=- 1 2.  To  do  this  we  will  substitute  - 1 2  for  A,  in  |  A' 

and  solve  the  matrix  equation  A'  X  =  0  ,  When  we  do  this,  we  get 
the  following  augmented  matrix. 


( 1 

0 

0 

0 

0 

V3 

1 

0  ' 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

-VI 

1 

0 

0 

0 

0 

1 

0 

- 1 

1 

0 

0  ' 

0 

0 

0 

1 

2 

1 

0 

u 

0 

0 

0 

0 

0 

1 

0  J 

This  augmented  matrix  above,  can  be  interpreted  as  the  following 
equations. 
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If  we  let  y3=l.  then  the  resulting  eigenvector  is 


X 


X=-I2 


'  -  VI  ' 
0 

VI 

i 

-2 

.  1  > 


The  graphical  interpretation  of 


X 


X=-I2 


can  be  seen  in  Figure  3.10. 


System  Motion: 


Figure  3.10 


orthogonal  set  of  eigenvectors  of  A  associated  with  the  eigenvalues 
0,0, 0,-6, -6  and  -12,  respectively.  We  normalize  these  orthogonal 
vectors  to  get  the  following  orthonormal  set  of  eigenvectors. 
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If  we  let  the  eigenvectors  above  form  the  columns  of  a  matrix  P, 

- 1  -Ik  k 

then  P  is  invertible  and  P  AP=D,  so  P  - — AP  =  — D  where  D  is  a 

4m  4m 

diagonal  matrix  with  the  eigenvalues  0,  0,  0,  -6.  -6,  and  -12  as  the 
entries  on  the  diagonal.  Since  our  goal  is  to  solve  the  differential 

equation  X  =  A  X  .  vve  will  let  U  =  P  X  ,  then  apply  Exercise 

2.4  with  n=6,  where  we  have  factored  — out  of  the  matrix  A 

4m 

„  k  _  k 

instead  of  — .  Thus  P  U  =  -;^P  D  U  becomes 
m  4m 


n(l)  -  n(2)  . 

P  Uj  +  P  0^  + 


..  +P 


P  X,U,  +  P  X^u^-H 


. .  +  P  u 
6 


where  Xj  =X2=X2=0,  X^=X5  =  -6.  and  Xj  =- 1  2.  This  equation  can  be 

simplified  by  multiplying  through  by  ,  gathering  terms  and 
moving  everything  to  the  left  side. 
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D*'*.-.  k  n<l\  1  rn>2)-  k  r.®,  1 

r  “r+srP  ^“iJt  “2-^^  ^2“2j-^ 


Jp‘«Q  „  1=^ 

L  6  4rn  6  6j 


Factoring  out  ^  \  . P^^\  we  have 


Since  the  columns  of  P  are  orthonormal  eigenvectors  of  A,  we  know 
they  are  linearly  independent.  Thus,  we  have  a  finite  linear 
combination  of  linearly  independent  vectors  which  equals  zero,  so 
the  coefficients  of  P  \  P^^^ . P^^^  must  be  zero.  If  we  set  each  of 

the  coefficients  in  the  above  equation  equal  to  zero,  we  have 


u  u  =  0 

I  4m  1  I 


u  =  u 

1  4m  I 


“2-4srS^=‘> 


or  u.  =  ^X_u. 

2  4m  2  2 


u  .  -  —  X  .u  .=  0 
6  4m  0  6 


%=4T^V6 


Since  Xj  =X9=X2=0,  =  and  X^=- 1  2.  these  differential  equations 


become 
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80 


X  =P  u 


V  J  ri  V  J 

^  /r 


3 
3 
3 


0  Y 

6 

1 

2 


'/I  :Zil 

6  7 

73 


0  —  --  - 

V/  ^  ^ 


0 


0 

3 

3 

3 


JK 

14 

/i4 

7 

14 

Th  774 


0 

Til 

7 


0 


0  —  -  - 


14 

Jtt 

14 


_J_ 

2 

0 


1 

2 

Jl 

6 

3 

714  73 

7  6  ; 


0 


C,,t-KC 


12 


Si' -^^22 

Si'-^S2 


c^,  cos 
41 


C51COS 


(^0 

(7¥‘) 

(7¥0 


+  c^_sin 
42 


+  C52S.n 


+  c  ._sin 
62 


(7¥  0 
(71 0 

(7¥  0 


The  solution  X  to  X  =-: —  A  X  is  found  by  multiplying  the  matrix  P 

4m 

by  the  vector  U.  The  components  of  the  solution  vector  X  are 


x.  =  c..^t  +  c..^-c..7^t-c,  ^ 


I  113 


12  3  31  6 


'32  6 


VTT_. 


VIT 


+  c  ,2^sm  t)-  c„.-:-^cos 


51  14 


(71  •) 


-  C 


52 


(71-  01  (7¥  0- =62  i  (7¥  ■) 


Y  =c  ^f4r  ^  +  r  ^ 

2  II  3  12  3  31  3 


t  +C 


VI 


+  C 

32  3  51 


#™(71') 


-  C 


5l 


I 

I 

i 
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It  is  possible  to  determine  the  values  of  the  Cjj  provided  we  have 

been  given  enough  details  about  the  system  which  we  have 
modeled.  Since  the  values  of  m  (the  mass  of  the  spring)  and  L  (the 
length  of  the  spring  in  equilibrium)  are  given,  and  if  we  know  the 
value  of  each  Cjj,  then  we  will  be  able  to  find  the  value  for  each  Xj 

and  Vj  at  a  given  time. 

We  now  want  to  apply  the  theory  of  the  spring-mass  system 
which  we  have  just  studied  to  understand  how  this  system  models 
the  vibrations  of  a  water  or  H2O  molecule,  as  shown  in  Figure  3.1  1 


Oxygen 


Figure  3.11 

Figure  3.1  has  three  springs,  but  the  water  molecule  has  only  two 
bonds.  The  third  spring  represents  the  repulsion  force  of  the  two 
hydrogen  atoms.  A  water  molecule  which  lies  in  the  xy-plane  would 
have  a  translational  motion  in  both  the  x-  and  y-directions  as 
described  in  Figures  3.5  and  3.6.  Also,  the  molecule  would  be  able 
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to  rotate,  as  we  saw  in  Figure  3.7.  Moncrief  and  Jones  (1977) 
explain  the  three  vibrational  modes  for  H2O  using  Figure  3.12. 
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The  vibration  in  Figure  3.12(a)  is  called  a  symmetric  stretch  since 
the  bonds  between  the  both  hydrogen  atoms  and  the  oxygen  atom 
are  stretching  by  the  same  amounts  at  the  same  time.  We  have 
already  considered  this  type  of  motion  in  Figure  3.8.  In  Figure 
3.12(b)  we  see  an  asymmetric  stretch  which  is  due  to  the  fact  the 
bonds  between  the  hydrogen  atoms  and  the  oxygen  atom  are  being 
stretched  by  the  same  amount  but  not  at  the  same  time.  Figure  3.9 
describes  th’is  same  mode  of  vibration  in  the  closed  spring -mass 
system.  The  final  mode  of  vibration  is  called  symmetric  bending 
and  is  seen  in  Figure  3.1  2(c)  This  vibrational  mode  consists  of  the 
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hydrogen-oxygen  bonds  remaining  at  the  same  length,  but  the  two 
hydrogen  atoms  vibrate  by  moving  further  apart  then  closer 
together.  We  have  already  seen  this  in  Figure  3.10. 
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Appendix  A:  Review  of  Differential  Equations 

We  will  limit  our  discussion  to  second  order  linear  differential 
equations  with  constant  coefficients.  This  appendix  is  not  meant  to 
replace  a  differential  equation  course,  but  only  to  show  how  to  solve 
a  very  select  group  of  differential  equations.  The  second  order 
linear  differential  equations  which  we  want  to  solve  are  of  the  form 

(A.i )  x+  m  x=  0. 

where  the  coefficient  of  the  x-term  is  a  constant  which  we  denote 
by  m.  Any  second  order  differential  equation  which  can  be  put  in 
the  form  of  Equation  (A.I)  is  called  a  linear  differential 
equation.  The  differential  equation  x  +  m  sinx  =  0  is  no  longer 
linear  because  sinx  is  a  nonlinear  function  of  x.  The  method  used  to 
solve  the  differential  equation  (A.  1 )  above,  depends  on  the  value  of 
m.  We  will  consider  three  possible  cases. 

CASE  1.  m=0 

If  m=0,  then  our  second  order  differential  equation  becomes 

x=  0. 

By  the  Fundamental  Theorem  of  Calculus,  if  x=0,then  x  =  a.and 
x  =  at  +  b. 
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Conversely,  if  x  =  at  +  b .  then  differentiating  this  equation 

with  respect  to  time  we  have 

^  =  a  or  X  =  a 
dt 

Differentiating  again,  we  obtain 

■^  =  0  or  x=0 
dt 

Therefore,  we  conclude  that  x  =  at  +  b  is  the  solution  to  the 
differential  equation  x=0. 

CASE  2.  m<0 

The  differential  equation  x  +  m  x  =  0  can  be  rearranged  as 
X  =  -  m  X  where  -m  is  a  positive  number.  Recall  from  calculus,  that 
the  exponential  function,  when  differentiated,  yields  a  multiple  of 
itself.  Thus,  we  want  an  exponential  function  which  when 
differentiated  twice  results  in  a  positive  multiple  of  itself.  Let  us 
pause  for  a  moment  and  consider  two  examples  of  exponential 
functions. 


x  =  and  x  =  e~^^ 

Taking  the  first  derivative  of  these  two  function  with  respect  to 
time,  we  obtain 
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x=2e^^  and  x  =  -2e"^^. 

After  taking  the  second  derivative,  we  have  the  following  two 
functions 

x=4e^^  and  x  =  4e“^V 

If  we  substitute  the  values  for  x  and  x  into  the  differential  equation 
X  -  4x  =  0 ,  we  see  that  x  =  e^^  and  x  =  e~^*  are  both  solutions  to  the 
same  differential  equation.  Furthermore,  any  linear  combination  of 

2t  -2t  . 

these  two  solutions  such  as  x=  CjC  ^-c^e  ,  is  also  a  solution  to 

X  +  mx  =  0  when  m=-4.  From  this  we  conclude  that 

y-m t  -V-m t  .  .  „ 

xsc^e  ■‘■^2^  ,  IS  a  solution,  for  all  Cj  and  C2. 


CASE  3.  m>0 

The  differential  equation  x+  mx  =  0  can  be  rearranged  as 
X  =  -  rn  X  where  -m  is  a  negative  number.  Recall  from  calculus,  that 
the  cosine  function,  when  differentiated  twice  yields  a  negative 
multiple  of  itself.  This  is  also  true  for  the  sine  function.  Let  us 
pause  for  a  moment  and  consider  two  examples  involving  the  cosine 
and  sine  functions. 


x=cos2x  and  x  =  sin2x 
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Taking  the  first  derivative  of  these  two  function  with  respect  to 
time,  we  obtain 


x  =  -2sin2x  and  x=2cos2x. 

After  taking  the  second  derivative,  we  have  the  following  two 
functions 

x  =  -4cos2x  and  x  =  -4sin2x. 

If  we  substitute  the  values  for  x  and  x  into  the  differential  equation 
X  +  4x  =  0 ,  we  see  that  x  =  cos2  x  and  x  =  sin 2x  are  both  solutions  to 
the  same  differential  equation.  Furthermore,  any  linear  combination 
of  these  two  solutions  such  as  x  =  CjCos2x  +  c  ^sin2x ,  is  also  a 

solution  to  X  +  mx  =  0  when  m=4.  From  this  we  conclude  that 
x  =  CjCos  Vrir t+  c^sinVirT  t,  isa  solution,  for  all  Cj  and  C2- 

Just  as  it  was  shown  in  Case  I,  where  m=0.  every  solution  of  x  =  0 
must  be  in  the  form  at*b.  It  can  also  be  shown  that  every  solution 
for  Cases  2  and  3,  where  m-0.  must  be  in  the  forms  we  have 
presented. 
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Appendix  B:  Solutions  to  Exercises 


=  D  U  can  be 


Exercise  2.1 

The  system  U  = 

rewritten  as  the  following  two  second  order  differential  equations 

-3k 


r-,1 

-  ^  1 

^  -  3 

f"'l 

m  1 

<  0 

-1  J 

m 


(B.l) 


u^  = 


^u 
m  '^2 


Similarly,  the  system  X  = 


f''  1 

=-( 

-2 

n 

[’‘'1 

m  \ 

.  1 

-2) 

U2J 

=  ^  A  X  can 


be  rewritten  as  the  following  two  second  order  differential 
equations 


(B.2) 


_  _ 2k  k 

^r-rTT^  +  nT  >^2 


Y  -  A  Y  2k 

2~  m  i  m  ^2 


Each  equation  in  (B.l)  can  be  solved  independently  using  only  basic 
techniques  from  differential  equations.  However,  since  each 
equation  in  (B.2)  is  in  terms  of  both  variables  Xj  and  X2<  neither 


equation  can  be  solved  independently.  Thus,  it  is  much  easier  to 
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solve  the  system  of  differential  equations  given  by 


1  L  k  - 


uj  m 


f-3  O'!  ^ 

I  0  -  I  J  u, 


r 


We  have  already  found  the  eigenvalues  of  A  to  be  Xj  =-3  and  X2=- 1 , 
3k  Ik 

thus,  r,  = —  and  ,  which  are  both  greater  than  zero.  From 

1  m  ^  m  ^ 

this  we  see  that  the  two  second  order  differential  equations  in  Uj 

and  U2  are 


-# 


These  are  both  second  order  linear  differential  equations  which  can 
be  solved  using  basic  techniques.  Their  solutions  are 


u  I  =  c  1 1  cos 


U2=  c^iCos 


(7¥  0*s2^'"^^/¥0 
(VI’  0*‘=22^'"(V¥0 


The  solution  to  the  original  system  of  differential  equations 
^  k 

X  A  X  is  found  by  substituting  the  values  for  both  the  matrix  P 
and  the  vector  U  into  the  equation  X  =  P  U  . 


f 


I 


I 
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Multiplying  the  matrices  on  the  right  side  together  and  equating 
components,  we  get  the  following  solutions  to  the  differential 
equation  that  models  the  spring-mass  system. 


Exercise  2.3 

Since  we  have  n  variables  Xj ,  X2 . we  have  n  equations  of 

motion  which  are 


d 

r  aT  1 

av 

d 

r  aT  1 

av 

d 

■  aT  1 

av 

dt 

ax,’ 

dt 

■;  ’  dt 

3Xn 

The  easiest  way  to  construct  these  equations  is  to  find  each 
component.  To  find  the  left  side  of  each  of  the  equations  of  motion, 
we  first  differentiate  the  equation  for  kinetic  energy  with  respect  to 
X.  (i«  1 , 2 . n). 
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+  0  +  ...+  Oj=  mx  I 


||-  =  im[0  +  2X2+0  +  ...+ 


-|^  =  jm[0  +  ...+  0  +  2x„]=mXn 
When  we  differentiate  each  of  these  with  respect  to  time,  we  have 


The  right  side  of  each  of  the  equations  of  motion  is 
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^  =  -ik[2b„x,+  2b,2X2  +  ...+  2b,„x„]  =  -k[b,,x,+b,2X2  +  ...-^b,_,x„] 


•^  =  4k[2b,2X,  +  2b22’‘2-"--"2b2„x„]=-k[b,2X,+b22=‘2-"--^'’2n’‘n] 


av  i 


=  in* ,  *  2b  2„X2+  ...+  2b„„x„]=-k[b,_,x,  +b2„X2+  ...  +  b„„x„] 


If  we  combine  these  components,  the  equations  of  motion  become 

^  =![-*>  n^-b, 2*2- 

k 


:  =  rn[-b,2X,-'’22X2-" 


b ,,  xl 
2n  nj 


>'n  =  rn[-'’-  x  -'’-  X--- 


In  1  2n  2 


-b„„X  1 
nn  nJ 


The  equations  of  motion  can  be  rewritten  in  matrix  form  as 


X  = 


f  X  ^ 

f-b  — b  ...  -b  ^ 

1 

II  12  In 

^2 

_  k 

-b,2  -<y^  ...  -b^^, 

-  m 

: 

k  > 

^  "^In  ~^2n  -  "‘^nn  , 

V 

m 


B  X 


95 


where  B  is  an  n  x  n  symmetric  matrix.  Now  we  have  an  equation 
that  should  look  very  familiar  to  us. 


Since  B  is  a  symmetric  n  x  n  matrix,  there  exists  an  orthogonal 
matrix  P  such  that  P  *BP=D.  The  matrix  D  is  the  diagonal  matrix 

whose  entries  along  the  main  diagonal  consist  of  the  eigenvalues 
Xj . Xf,  of  B  and  the  columns  of  P  are  the  corresponding 

— >  — > 

eigenvectors  X  . X  associated  with  the  eigenvalues  X  ...  .  Xp , 

respectively.  The  following  theory  will  be  very  similar  to  the 
theory  that  we  developed  for  the  spring-mass  system  with  two 
blocks,  except  the  sizes  of  the  matrices  and  vectors  will  be  n  x  n  and 
n  X  1,  respectively.  If  we  multiply  both  sides  of  this  equation  by  P  * 

and  use  the  identity  P  P  ^  =  Ip,,  we  get 


(p  P  0  X  =  (p  ' 


_i 

To  simplify  this  equation,  we  let  U  =  P  X.  To  introduce  the  vector 
variable  U ,  we  substitute  U=P  X  into  this  equation.  We  then 

- 1  -»  _i 

substitute  P  BP  =  D  and  U  =  P  X  into  the  right  side  to  obtain 
^  k 

U  =  -pfr  D  U  .  Now.  multiplying  both  sides  of  this  matrix  equation  by 
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the  matrix  P.  we  get  P  U  =  ^P  D  U .  Rewriting  the  left  side  of 


P  U  =-^P  D  U  gives 


r.(2)  „(n)^  U  (I).  (2).  „(n)  _ 

PU=VP  P  ...P  )  2  -p  u,  +  P  u  +...  +  P  u 

:  I  ^ 


where  P  represents  the  ith  column  (i  =  1.2 . n)  of  the  matrix  P. 

We  rewrite  the  right  side  to  obtain 


irPD  P“>-n 


Thus,  P  U  =  ^P  D  U  can  be  written  as 


Jl)  -  „(2)  .  ^(n)  ..  k  fJl),  „(2),  ^(n),  1 

P  u,  +  P  u^+.-.+  P  Un=7^rLP  X^u.^+...+  P  X^UnJ 


We  can  simplify  this  equation  by  multiplying  through  by  ~  , 
gathering  terms,  moving  everything  to  the  left  side,  and  factoring 

r,(l)  r>(2)  ^(n)^ 

out  P  ,  P  . P  from  each  quantity. 


[“,  -  m  ^1 'J|]p"’+[“2' ■^^2  “2] ^  [“n  -  ‘'nip'"’ =  0 
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Since  the  columns  of  P  are  orthogonal,  we  know  they  are  linearly 
independent,  thus,  the  coefficients  of  the  vectors  \  ...  .  P* 


must  be  zero.  If  we  set  each  of  the  coefficients  in  the  above 
equation  equal  to  zero,  we  have 


-nr^n“n 


These  are  all  second  order  linear  differential  equations  which  can  be 

k. 

solved  using  basic  techniques.  If  we  let  f'j  =  ~  TfT^j .  where  i=  1 ,  2 . 

n.  then  these  equations  become 

Uj+r|Uj  =  0.  U2  +  r2U2  =  0 . Un  +  rnUn=0 


Using  the  following  formulas,  we  can  solve  for  the  vector  U.  (Note: 
To  determine  whether  rj  is  zero,  negative  or  positive,  substitute  Xj 

into  r.  =- -^X. .) 

1  m  1 


If  r.  =  0,  then  u.  =  c.,t  +  c  „ 

I  I  tl  i2 


r  t  -r  .i 

If  r.<0,  then  u.  =  c,.e  '  +c  „e  ' 

I  I  il  i2 
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The  solution  to  the  original  system  of  differential  equations 
^  k 

X  =-^  B  X  is  found  by  substituting  the  values  for  both  the  matrix  P 
and  the  vector  U  into  the  equation  X  =  P  U . 

Exercise  ZA 

Since  the  spring-mass  system  lies  free  in  the  xy-plane.  the  entire 
system  can  move  vertically  up  or  down,  horizontally  to  the  left  or 
right,  or  rotate.  These  types  of  motion  are  called  rigid  motions.  The 
system  can  also  vibrate  producing  the  motions  that  are  described  by 
Figures  2.6  and  2.7. 

Exercise  2.5 

(a)  Suppose  the  three  masses  are  moved  to  the  right  causing  the 
first  three  springs  to  stretch  by  different  amounts  and  causing  the 
fourth  spring  to  be  compressed.  This  is  depicted  in  Figure  B.l. 


m  m  m 


l-x,H  hXz-H  l-XjH 


Figure  B.  I 

First,  we  need  the  equation  which  describes  the  kinetic  energy  of 
the  system  in  Figure  B.  I . 


99 


^  1  f.  2  .  2  .  21 

The  potential  energy  of  the  system  is  the  sum  of  the  potential 
energies  of  each  spring.  Spring  1  is  stretched  from  its  equilibrium 
position  by  the  amount  Xj.  so  the  potential  energy  for  spring  1  is 

V  j  X  j .  Spring  2  is  stretched  from  its  equilibrium  position  by 

the  amount  Xj-X2.  so  that  ”^2)  the  potential  energy 

for  spring  2.  Spring  3  is  stretched  from  its  equilibrium  position  by 
the  amount  X2-X3,  producing  a  potential  energy  of 

^3“  1  *^(^2*^3)  foJ' spring  3.  Spring  4  is  compressed  from  its 
equilibrium  position  by  the  amount  x^.  Thus,  the  potential  energy 

for  spring  4  is  V^  =  -7kx^.  Therefore,  the  potential  energy  of  the 

4  ^  o 


system  is 
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Since  we  have  three  variables  X| .  X2.  and  x^,  we  have  three 
equations  of  motion  which  are 


d 

■  dT  ■ 

BV 

d 

r  1 

dV 

anri 

r  aT  1 

av 

dt 

dt 

3X3 

dt 

^^3 

The  easiest  way  to  construct  these  equations  is  to  find  each 
component.  To  find  the  left  side  of  each  of  the  equations  of  motion, 
we  first  differentiate  the  equation  for  kinetic  energy  with  respect  to 
X.  (i=1.2.3). 


aT 

=  m  f- 

2X.  +  0  +  0' 

1=  m  X. 

3x, 

2  L 

1 

J  1 

3T 

=  4-  m  f 

0  +  2  X  .  +  0 

•1  =  m  X  _ 

ax^ 

2  L 

2 

J  2 

3T 

1  r 

0  +  0  +  2  X  ^ 

3X3 

=  ^m[ 

J  =  mx3 

Now,  differentiating  each  of  these  with  respect  to  time,  we  have 
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dt 
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ax 


U 


dt 


_dj_ 

ax. 


=  •^["’*2]=  ">*2 


dt 


ar 

ax. 


The  right  side  of  each  of  the  equations  of  motion  is 


•|^  =  -k[2x,-X2]=k[-2x,  +  X2] 


^  =  -k[2x2-x,-X3].k[x,-2X2H-X3] 

2 


^=-k[2X3-X3]=k[X3-2X3] 


If  we  equate  these  components,  th:  equations  of  motion  become 


*3  =  ■^[’‘2 -2’' 3] 
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The  equations  of  motion  can  be  rewritten  in  matrix  form,  which  is 
the  system  of  differential  equations  modeling  the  spring-mass 
system  in  Figure  B.l . 


X  = 


X.  ^ 


‘3  } 


JL 

m 


( -2 
I 

0 


-2  I 
I  -2 


X  > 
^1 

2 

/ 

X., 

^  3  J 

(b)  Since  A  is  a  symmetric  matrix,  all  of  its  eigenvalues  are  real  and 
A  is  diagonalizable.  We  begin  by  finding  the  eigenvalues  of  the 
matrix  A. 


f -2-X 


det(A  -XI)  =  det 

V 


1 

0 


1 

-2-X 

1 


=  (X  +  2)(X  +  2-V2)(X  +  2  + V2) 


XJ 


If  we  set  det(A  -XI)  equal  to  zero  and  solve  for  X.  we  find  the 
eigenvalues  are  Xj  =-2,  X^  =  -2*V~^ and  X2=-2-4~2.  Thus,  there  exists 

an  invertible  matrix  P  such  that  P'*  AP  =  D.  D  is  the  diagonal  matrix 
whose  entries  along  the  main  diagonal  consist  of  the  eigenvalues  of 
A  and  the  columns  of  P  are  the  corresponding  eigenvectors.  To  find 
P.  we  need  to  find  the  eigenvectors  associated  with  Xj  =-2,  X^  =  -2*V~^ 

and  X^’‘-2-'i~2.  For  Xj  =-2  we  have  to  reduce  the  following 


augmented  matrix 
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f -2+2  1  0 

1  -2+2  1 
0  I  -2+2 


0  ^ 
0 

0  y 


to  obtain 


^  1  0  1  I  O' 
0  I  0  I  0 

.  0  0  0  I  0  > 


which  yields  the  following  equations:  x,=- 


and  X2=0.  If  we  let  x  j  =- 1 ,  an  eigenvector  associated  with  Xj  =-2 


is  X 


^  -  1  ^ 
0 
I  J 


For  X2“"2+V~2* we  have  to  reduce  the  following  augmented 


matrix. 


^-2-(-2  +  V2)  1 

I  -2-(-2  +  V2) 


0 


0  I  0 

I  I  0 

-2-(-2  +  V2)  I  0 


10-10 


,  which  yields  the  following  equations: 


to  obtain  I  0  -  I  VJ  |  0 

0  0  0  I  0  J 

X  J  =X2  and  X2=^fYxy  If  we  let  x^=  I ,  an  eigenvector  associated  with 


Xj=-2>V2  is  X 


X=-2+/2 


f  1  ^ 

VI 


V  I  J 
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For  X,2=-2-V^ we  have  to  reduce  the  following  augmented  matrix. 

'-2-(-2-V2)  1  0  I  o' 

1  -2-(-2-V2)  1  I  0 

.  0  1  -2-(-2-V2)  I  0, 

"10  -  1  I  0  " 

to  obtain  o  -1  -VJ  |  0  ,  which  yields  the  following 

.  0  0  0  I  0  > 

equations:  x  j  =X2  and  ^3"  •  •  eigenvector 


associated  with  A.,  =-2-V2  is  X  _=  -^/2  ■  Using  the  theory 

*  X  =-2  -J2 

3  V  1  y 

that  we  developed  in  Exercise  2.3,  we  know  that  we  must  first  solve 

k 

the  differential  equation  P  U  =  -^P  D  U .  which  leads  to 


K  -  w  s «,]  p'  “  2  - 1^2  “2]  p"’  ^  3  - “3]  p™  =  ° 


Since  the  columns  of  P  are  orthogonal,  we  know  they  are  linearly 

independent.  Thus,  the  coefficients  of  the  vectors  P\s((  1 ),  ).  P\s((2), 
(3) 

)  and  P  ,  must  be  zero.  If  we  set  each  of  the  coefficients  in  the 
above  equation  equal  to  zero,  we  have 
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Substituting  in  the  eigenvalues,  we  obtain 

02-|-(-2  +  V2)u2  =  0.  u3-|-(-2-V2)u3 


After  simplifying  these  equations  they  become 


+  =0'  ^2* 


(2-V2)k 


U2=0, 


^3^ 


(2  + VDk 


From  these  differential  equations  we  observe  that 

_2k  (2-V2)k  (2  V2  )k 

*'l  ~  m  ’  ^2~  ’  ^3~  m 

Since  each  rj  (i=l,  2,  3)  is  greater  than  zero,  the  solutions  to  these 

second  order  linear  differential  equations  are 


u,=C||Cos( 

(  /(2- V2)k  .V  . 

1 

CM 

“2  ==2,  cos 

m  tj  +  C22Sin 

W  m 

(  /(2  + V2)k  ^ 

f  /(2  +  V2)k  ^ 

“3=  =3,005 

W  m  V'^^32^"^ 

K^l  m  U 
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The  solution  to  the  original  system  of  differential  equations 
^  k 

X  =-^  A  X  is  found  by  substituting  the  values  for  both  the  matrix  P 
and  the  vector  U  into  the  equation  X  =  P  U  . 


Multiplying  the  matrices  on  the  right  side  together  and  equating 
components,  we  get  the  following  solutions  to  the  differential 
equation  which  model  the  spring-mass  system. 
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X2=‘n/2  C^jCOS 


-V2  C3,cos(7^^*1^—  t)-V2  CjjSin  t) 


X3  =  c,|Cos(y^  t)  +  C|3Sin(y^t) 


+  C2jCOsl 

u 

f(2-V2)k  ' 

Uc22sinl 

/(2-V2)k  ) 

m  V 

Iv  m  ‘•J 

+  C3,cos 

f(2  +  V2)k  .1 

m  '■J 

kc32sin 

(c)  To  describe  the  configurations  in  which  the  spring-mass  system 
vibrates,  we  need  to  write  each  Uj  in  terms  of  the  Xj.  This  can  be 

— ♦  — >  — >  _j  — > 

done  by  using  the  matrix  equation  X=PUor  U  =  P  X.  We  first 
find  P  *  and  substitute  it  into  the  matrix  equation. 


(  1 

0 

1  ^ 

1 

2 

2 

r  ^  1 

1 

u 

1 

V2 

1 

1 

X 

2 

4 

4 

4 

2 

V  ^3  > 

1 

V2 

1 

^  ^3  > 

1  4 

4 

4  ; 

Multiplying  these  matrices  together  and  equating  components,  we 
get  the  following  three  equations. 
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2 


V2 


X,  +  V2  X2  +  X2 


*^2= - 4 


X  - V2  X  +x 
^3“  4 


These  three  equations  describe  the  relationships  between  the 
variables 

Xj.  X2  and  x^.  The  three  modes  in  which  this  spring-mass  system 


-X  +x 

vibrates  are  u^  = - — — -  with  U2=*0  and  u^'O. 


X  +V2  X  +x 

u^=  — - r — - -  with  u,  =0  and  u,=0.  and 

2  4  I  j 


X  -  VJ  X  +x 

u^=  — - - - — - - -  with  u 

3  4 


=0  and  U2=0. 


In  the  first  mode.  U| 


-  Xj  +  x 


3 


2 


with  U2=0  and  u^^O.  Since 


-X  +x 

- - — -  represents  how  the  distance  between  blocks  one  and 


three  is  changing,  the  first  mode  of  vibration  describes  how  the 
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distance  between  the  two  blocks  is  changing.  We  visualize  this  by 
considering  a  series  of  diagrams  similar  to  those  in  Figure  2.6. 
Recall,  the  banner  is  made  of  an  elastic  material  and  indicates  the 
distance  between  the  two  blocks.  When  this  system  vibrates,  we 
see  the  banner  contracting  and  stretching  with  a  frequency 
associated  with  Xj.  This  is  indicated  by  the  series  of  “snapshots"  of 

the  spring-mass  system  in  motion  seen  in  Figure  B.2. 


Figure  B.2 


In  the  second  mode, 


1 10 


X  +  V2  X  +x 

— = - - — with  u ,  =0  and 

4  1 


0^=0.  To  visualize  this  mode  of  vibration,  we  consider  a  series  of 

diagrams  similar  to  those  in  Figure  2.6.  When  this  system  vibrates, 
we  see  the  flag  moving  back  and  forth  with  a  frequency  associated 
with  X,2.  This  is  indicated  by  the  series  of  “snapshots”  of  the  spring- 

mass  system  in  motion  as  seen  in  Figure  B.3. 


Figure  B.3 


1 1 1 


In  the  third  mode.  u^  = 


X  -  VJ  X  +x 

— - : — with  u ,  =0  and  Uo=0. 

4  I  ^ 


To  visualize  this  mode  of  vibration,  we  consider  a  series  of  diagrams 
similar  to  those  in  Figure  2.6.  When  this  system  vibrates,  we  see 
the  flag  moving  back  and  forth  with  a  frequency  associated  with 

This  is  indicated  by  the  series  of  “snapshots"  of  the  spring-mass 
system  in  motion,  as  seen  in  Figure  B.4. 


1  12 


Exercise  2  6 

Since  the  spring-mass  system  lies  in  the  xy-plane.  the  entire  system 
can  move  vertically  up  or  down,  horizontally  to  the  left  or  right,  or 
rotate.  Besides  these  rigid  motions,  the  system  can  also  vibrate 
producing  the  motion  that  were  described  in  Exercise  2.6. 


Exercise  3  1 

We  write  f  as  a  Taylor  series  expansion  expanded  about  zero. 


(^)=r(o) 


r  (0)  I  r  "(o)  I'l 

M  L 


We  use  f  which  has  been  repeated  below  for  convenience,  to 
find  f(0).  r(0)  and  r  (0). 


f 


-X2)-V3(y,-yJ)+i((, 


f(0)  =  1 

r(0)  =  l[(x|-X2)-V3(y|-y2)] 


r  ■(0)  =  -7[(x,-X2)-(y,-y2)P+(x,-X2)  -(y|-yj 


Substituting  these  into  the  Taylor  series  expansion  and  simplifying, 
we  have 

IlJ"  *  I!  L 


-xJ-V3(y,-y2)f 


+  higher  power  terms. 


Exercise  3.2 

Figure  B.5  contains  the  coordinates  of  the  blocks  and  the  distance 
formula  for  the  spring  from  Figure  3.3. 


Figure  B.5 
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The  potential  energy  of  this  spring  is  one  half  the  spring  constant  k 
times  the  square  of  the  distance  that  the  spring  is  stretched.  The 

displacement  of  the  spring  from  its  equilibrium  position  (the 
distance  that  the  spring  is  stretched)  is  |d-L|.  Expressing  the 
potential  energy  of  the  spring  in  terms  of  |d-L|.  we  have 

V23=ik|d-L|2 


We  want  to  rewrite  using  the  variables  X2.  ^2'  ^3  y3- 


do  this,  we  must  first  simplify  the  expression  for  the  distance  d. 

1 


2 


1 


We  can  rewrite  the  quantity  expressed  by  the  square  root  using  its 
Taylor  series  expansion.  The  terms  of  higher  powers  have  been 
grouped  together  for  convenience. 
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d  =  L 


1  + 


1 ! 


+  terms  of  higher  powers 


This  expression  can  be  simplified  by  multiplying  through  by  L  and 
then  moving  the  L  to  the  left  side.  The  resulting  quantity  is  what 
we  wanted  to  find. 


d-L  =  j  (^2  higher  powers 


This  quantity  can  now  be  substituted  into  the  formula  for  the 
potential  energy  V23. 


V33  =  ik|d-L 


+  terms  of  higher  powers 


1 


The  Taylor  series  expansion  for  potential  energy  can  not  have 
nonzero  linear  terms,  since  the  system  has  an  equilibrium 
configuration.  Also,  we  are  considering  only  smaM  vibrations  so  we 
ignore  the  terms  of  higher  powers.  The  formula  for  potential 
energy  V23.  is 
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''23=214  [(>‘2-’‘3)+^(y2-y3)] 


k. 

2 


22^2.2  ^  ^ 
4'*’4  4  4  2  ^2^3^  2  ^2^2  2 


^^3^2 


Vb  v^  _^v 

2  2^3  2  3^3  2^2^3 


J 


Exercise  3.3 

Figure  B.6  contains  the  coordinates  of  the  blocks  and  the  distance 
formula  for  the  spring  from  Figure  3.4. 


(-^-3-  ^3) 

m-^ — yww 


X. 


where  d  = 


Figure  B.6 


The  potential  energy  of  this  spring  is  one  half  the  spring  constant  k. 
times  the  square  of  the  distance  that  the  spring  is  stretched.  The 

displacement  of  the  spring  from  its  equilibrium  position  (the 
distance  that  the  spring  is  stretched)  is  |d-L|.  Expressing  the 

potential  energy  of  the  spring  in  terms  of  |d-Lj,  we  have 


We  want  to  rewrite  V  j  ^  using  the  variables  x  j ,  y  j ,  and  y^.  To 


do  this  we  must  first  simplify  the  expression  for  the  distance  d. 


We  can  rewrite  the  quantity  expressed  by  the  square  root  using  its 
Taylor  series  expansion.  The  terms  of  higher  powers  have  been 
grouped  together  for  convenience. 


d  =  L 


I  + 


I  ! 


l|-+ higher  power  terms 


i 
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This  expression  can  be  simplified  by  multiplying  through  by  L  and 
then  moving  the  L  to  the  left  side.  The  resulting  quantity  is  what 
we  wanted  to  find. 

d  -  L  =  I  -  X2j+  higher  power  terms 


This  quantity  can  now  be  substituted  into  the  formula  for  the 
potential  energy  V  j  ^ 


V,3=ik|d-Ll' 


+  higher  power  terms 


The  Taylor  series  expansion  for  potential  energy  can  not  have  any 
nonzero  linear  terms,  since  the  system  has  an  equilibrium 
configuration.  Also,  we  are  considering  only  small  vibrations  so  we 
ignore  the  terms  of  higher  powers.  The  formula  for  potential 
energy  Vj 3,  is 


Exercise  3.4 
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Substituting  X 


X  =0 
I 


r  1  ^ 
1 
1 

0 
0 

vo; 


into  V 


-k 


X 

X  =0 

V  j  y 


=  X,  n-i^AX  =0. 

^  =0  8  m  X  =0 

J  j 


gives 


V 


'  *  1 

- 

1 

. 

1 

1 

=  (1110  0  0)  A 

1 

0 

8m 

0 

0 

0 

Lo> 

U  J 

^  -5 

1 

4 

VI 

-VI 

0  ' 

"  1  ^ 

1 

-2 

1 

-VI 

0 

VI 

1 

=  — -  (  1  1  1  0  0  o) 

4 

1 

-5 

0 

VI 

-VI 

1 

8m  ‘  ^  ^ 

VI 

-VI 

0 

-3 

3 

0 

0 

-VI 

0 

VI 

3 

-6 

3 

0 

^  0 

VI 

-VI 

0 

.  3 

-3  > 

1  0, 

-k 
8  m 


(  0  0  0  0  0  0) 


^  i 
1 
I 

0 
0 

0  J 


=  0 
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Exercise  3.5 


Substituting  X 


X  =0 
2 


^0^ 
0 
0 
1 
1 

V 1 ; 


into  V 


f-* 

X 


“^T  _  ic  -♦ 

=  X,  *^AX  =0, 
^  =0  8  m  X  =0 
J  j 


gives 


(0^ 

r  0  ^ 

0 

0 

0 

=  (  0  0  01  I  l)  A 

0 

1 

8m 

1 

1 

1 

^  1  > 

1 1  > 

8  m 


^  -5 

I 

4 

VI 

-VI 

0  ' 

"  0  " 

1 

-2 

1 

-VI 

0 

VI 

0 

1  1  1) 

4 

1 

-5 

0 

VI 

-VI 

0 

VI 

-VI 

0 

-3 

3 

0 

1 

-VI 

0 

VI 

3 

-6 

3 

1 

.  0 

VI 

-VI 

0 

.  .  3 

-3  > 

.  1  > 

-^(  0  0  0  0  0  0) 
8m 


^  0  > 
0 
0 

1 
I 

V  1  / 


=  0 


I2i 


-  5  1  4  VI  -VI  0  ' 

1  -2  1  -VI  0  VI 

4  1-5  0  VI  -VI 

VI  -VI  0-330 
-VI  0  VI  3-6  3 

0  VI  -VI  0  3  -3  > 


=  -^(-V^  ^VT-VX -3.0.3) 


T  x 


To  show  that  X  and  X  are  orthogonal,  we  show  their  dot 

>-,=-6 

product  is  zero. 


